DUDLEY U'-'TOT. LIEB.AB.Y 
HAA 7 AL P03'~'3 &A'0 TT . fE 
MOMTEPBY OAU*^- 



SCHOOL 

' 35345-BOOS 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

THE CURVE FITTING 
OF 

PORTABLE WEAR METAL ANALYSIS 
by 

Byung H. Min 
March 1987 

Thesis Advisor Harold J. Larson 



Approved for public release; distribution is unlimited. 



T233312 




UNCLASSIFIED 

SECURITY Class. FiCaNQn Qp This pag£ 



REPORT DOCUMENTATION PAGE 



la REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 



’b RESTRICTIVE MARKINGS 



2a SECURITY CLASSIFICATION AUTHORITY 



2b OEClASSiFiCATiON / DOWNGRADING SCHEDULE 



3 DISTRIBUTION/ AVAILABILITY OF REPORT 

Approved for public release; 
Distribution is unlimited 



4 PERFORMING ORGANIZATION REPORT NUM9ER(S) 



S MONITORING ORGANIZATION REPORT NUVBER(S) 



6a NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



6b OFFICE SYMBOL 
(If applicable) 

55 



7a NAME OF MONITORING ORGANIZATION 

Naval Postgraduate School 



Sc ADDRESS {C/fy. State, and ZlPCode) 

Monterey, CA 93943-5000 



7b ADDRESS (City State and ZIP Code) 

Monterey, CA 93943-5000 



8a NAME OF FUNDING /SPONSORING 
ORGANIZATION 



8b OFFICE SYMBOL 
(If applicable) 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



Be ADDRE SS (C/fy, State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 



program 


PROJECT 


TASK 


WORK UNIT 


ELEMENT NO 


NO 


NO 


ACCESSION NO 



'1 TITLE (Include Security ClitSifiCidon) 










THE CURVE FITTING OF PORTABLE 


WEAR 


METAL ANALYSIS 




12 PERSONAL AUTHOR(S) 

MIN, Byun^ 


Ho 










* 3 a TV PE OF REPORT 

Master's Thesis 


1 3b T’ME COVERED 
FROM _ 


TO 




14 DATE OF REPORT (Year, Month, Day) 

1987 March 


IS PAGE COoNT 

67 


’6 Supplementary notation 



17 COSATi CODES 


F ElD 


GROUP 


SUBGROUP 















18 SUBJECT TERMS ( Continue on reverie if neceisary and identify by block number) 

PWMA , Curve Fitting, Sum of Squares, Rational 
Equation, Spectroscopy. 



9 ABSTRACT (Conf*nu€ on revene if neceisary and identify by block number) 

This thesis describes the curve fitting algorithm that the Portable 
Wear Metal Analyzer used for calculating concentrations in ppm and 
compares this with some alternative algorithms. Each algorithm 
considered fits a curve to the three standards used, which were 20%, 
50%, and 100% of the full scale for all nine wear metals APL was 
used for all programs. 



20 0 S'R’BUT ON /AVAILABILITY OF ABSTRACT 

El UNCLASSIFIED/UNLIMITED □ SAME AS RPT □ DTlC USERS 


21 A8STRACT SECURITY CLASSIFICATION 

UNCLASSIFIED 


22a NAME OF RESPONSIBLE INDIVIDUAL 


22b TELEPHONE (Include Are* Code) 

408-646-247.3 


22 c OFFICE SYMBOL 

IxLd 



DO FORM 1473.84 MAR 83 APR edit.on may be u>ed until e*hau*ted SECURITY CLASSIFICATION OF Thi$ PAGE 

All other editions are obsolete 



i 



UNCLASSIFIED 



Approved for public release; distribution is unlimited. 



The Curve Fitting 
of 

Portable Wear Metal Analysis 



by 



Byung II. jVlin 
Major, Koredn Army 
B.S., Korea Military Academy, 1976 
B.S., Seoul National University, 1980 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OPERATIONS RESEARCH 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1987 



ABSTRACT 



This thesis describes the curve fitting algorithm that the Portable Wear Metal 
Analyzer used for calculating concentrations in ppm and compares this with some 
alternative algorithms. Each algorithm considered fits a curve to the three standards 
used, which were 20%, 50%, and 100% of the full scale for all nine wear metals. APL 
was used for all programs. 






TABLE OF CONTENTS 



-r/Tc f ^ 



L INTRODUCTION S 

II. PWMA AND DATA DESCRIPTION 9 

A. GENERAL 9 

B. THE PORTABLE WEAR METAL ANALYZER 10 

C. TEST PROGRAM AND DATA COLLECTED 11 

III. ALGORITHMS 16 

A. GENERAL 16 

B. ALGORITHM 1 17 

C. ALGORITHM 2 17 

D. ALGORITHM 3 20 

E. ALGORITHM 4 22 

F. ALGORITHM 5 23 

IV. ALGORITHM COMPARISON 27 

A. GENERAL 27 

B. SUM OF SQUARES METHOD 27 

V. CONCLUSIONS AND RECOMMENDATIONS 39 

APPENDIX A: PWMA FIELD TEST BURN SEQUENCE 40 

APPENDIX B: DATA PPM401 42 

APPENDIX C: PROGRAMMING LIST 43 

APPENDIX D: SUMMARY OF RANK 4S 

LIST OF REFERENCES 65 

INITIAL DISTRIBUTION LIST 66 



4 



LIST OF TABLES 



1. ABSORBANCE DATA OF TUBE 401 

2. DATA OF EACH LABORATORY 

3. TABLE OF OIL SAMPLE MIXTURES 

4. ESTIMATED CONCENTRATION BY 6 ALGORITHMS 

5. ALGORITHM COMPARISON FOR FE 

6. SUMMARY OF THE RANKS FOR TUBE 401 

7. RANKS OF 24 TUBES FOR FE 

S. RANKS OF 24 TUBES FOR AG 

9. RANKS OF 24 TUBES FOR AL 

10. RANKS OF 24 TUBES FOR CR 

11. RANKS OF 24 TUBES FOR CU 

12. RANKS OF 24 TUBES FOR MG 

13. RANKS OF 24 TUBES FOR NT 

14. RANKS OF 24 TUBES FOR SI 

15. RANKS OF 24 TUBES FOR TI 

16. THE BEST ALGORITHM FOR EACH ELEMENT 



12 

13 

15 

32 

32 

33 

34 

34 

35 

35 

36 

36 

37 

37 

3S 

38 



5 



mm mm m 



LIST OF FIGURES 



2.1 Calibration Technique used by the PW\1A 14 

3.1 An Asymptote Line in Algorithm 2 20 

.2 Possible Observed Curve from Algorithm 3 21 

.3 The Curve Bending Downward at Lower Concentration bv Algorithm 

4 7 '. . . 7 23 

.4 The Curve Bending Upward at High Concentration by Algorithm 4 24 

.5 The Curve Bending Downward at High Concentration bv Algorithm 

5 T '. . . T 25 

.6 The Curve Bending Upward at Lower Concentration by Algorithm 5 26 

4.1 The Horizontal Deviation of Residual by Algorithm 2 28 

4.2 The Horizontal Deviation of Residual by Algorithm 4 30 



6 



ACKNOWLEDGEMENTS 



In preparation of this thesis I am greatly indebted to Prof. Larson who gave me 
a great deal of valuable advice and 1 would like to express thank for his permission to 
use the data that I have employed. I found this was a good opportunity to have a good 
relationship between us. I am also indebted to Prof. Barr for his review of my draft. I 
want to express my gratitude to all my fellow friends who gave me good 
recommendations. I am gratefully thankful to my wife she always encouraged me and 
gave me a good confidence to do my job. 



7 



I. INTRODUCTION 



Spectrometric oil analysis is used to determine the type and amount of wear 
metals in lubricating oil samples. The Portable Wear Metal Analyzer (PWMA) was 
developed for the USAF to perform engine oil analysis when aircraft are deployed 
away from their home bases. The analyzer can detect nine wear metals in parts per 
million concentration ranges. PWMA takes samples with known concentration to 
obtain absorbance areas, and uses those to determine a calibration curve. Absorbance 
areas obtained from unknown samples are compared to the calibration curve so that 
their concentration can be calculated. 

This thesis describes the curve fitting algorithm that the PWMA used for 
calculating the concentration in ppm and compares this with some alternative 
algorithms for the curve fitting. We compare these algorithms to see which algorithm 
performs best. Each algorithm considered fits a curve to the three standards used 20%, 
50%. and 100% of the full scale for all nine wear metals. APL was used for all 
programs. 

Data gathered as the result of the field test of the PWMA are described in 
Chapter II. Three rational functions, and polynomials which are quadratic in 
absorbance and quadratic in concentration are suggested as alternatives in Chapter III. 
A squared error technique is employed to compare the algorithms' performances. 
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II. PVVMA AND DATA DESCRIPTION 



A. GENERAL 

The PWMA test program was conducted at Naval Air Station. Pensacola FL. 
and three Air Force field test sites, Langley AFB. VA, Myrtle Beach AFB, SC, and 
Elmendorf AFB. AK in 19S5. The purpose of this test plan was to determine the 
effectiveness and evaluate the functionality of the newly developed PWMA prototypes 
for measuring wear metal content in turbine engine lubricating oils under field 
conditions. 

When aircraft are deployed from their home bases to austere sites where normal 
oil analysis capabilities are nonexistent, there still exists a requirement to determine 
wear metal concentration in turbine engine lubricants. Due to the critical nature of 
turbine engine lubrication system component wear, particularly in high performance 
systems such as the F-15, F-16, and A- 10 aircraft, the requirement for oil analysis may 
be on an "after each flight" basis. Presently, the Air Force deploys Flame Atomic 
Absorption (FAA) spectrophotometers [Ref. 1: p. 50], which were not designed for 
mobile use. These instruments are sensitive and highly susceptible to damage when 
deployed for mobility exercises: such damage may lead to partial or complete loss of oil 
analysis capability in early stages of such exercises. In addition, the instruments require 
hazardous support accessories, such as compressed gases (nitrous oxide and acetylene) 
and solvents for dilution (and, in some past cases, acid solutions for digestion of wear 
metal particulates). 

In 19S0 F1Q TAC issued TAF SON 305-S0, establishing the requirement for a 
critically needed light-weight portable aircraft engine oil analysis unit. In response, 
AFWAL 1 undertook a contractual effort to determine the feasibility of developing such 
a device (Contract Number F33615-S0-C-2037). Since no such spectrometer with 
required light-weight portable features was commercially available, a major 
development program was initiated to ruggedize and miniaturize the conceptual 
hardware for mobility use. The actualization of the portable wear metal analyzer has 
been achieved under Contract Number F33615-S1-C-20S0 by Perkin-Elmer 
Corporation. Applied Science Division. 



l Air Force Wright Aeronautical Laboratory 
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B. THE PORTABLE WEAR METAL ANALYZER 

The PWMA features simultaneous measurement of nine wear metals in 
lubricating oils with an analysis cycle of four minutes (via air cooling): Silver (Ag), 
Aluminum (Al). Chromium (Cr). Copper (Cu). Iron (Fe), Magnesium (Mg), Nickel 
(Xi), Silicon (Si), and Titanium (Ti). Moreover, with the potential modification to 
water cooling on production units, this time can be reduced to two minutes. In 
addition, the PWMA may be operated on either 1 10 220 VAC, 50 or 60 Hertz, to 
accommodate use in both the L'.S and Europe. The PWMA is packed in two 
aluminum "suitcases": the electronics argon case, weighing forty pounds; and the 
furnace optics case, weighing sixty pounds. The cases are readily set-up and 
interconnected by two multiconductor cables and one pneumatic hose. Developed for 
simple operation and reduced dependence on skill and labor of the operator, the 
PWMA is characterized as follows: ease of set-up and self-contained configuration; 
minimal consumables (sample tips, graphite tubes and argon gas supply); no hazardous 
gases or chemicals (inert compressed gas only): microprocessor-controlled software for 
calibration, diagnostic etc.; and elimination of transcription errors by automated read- 
out, print-out of data and results. 

The PWMA operates on the principle of absorption of characteristic radiation 
by atoms of the wear metals of interest. Two multi-element hollow cathode lamps 
(HCLs), alternately pulsed for source combination and modulation, provide emission 
source lines for discrete line absorption by specific atoms, which are formed through 
heating in a preprogrammed ramped cycle (up to 3000° C) in the graphite furnace. All 
nine absorption lines can be separated simultaneously by the polychrometer and 
transmitted to respective discrete photomultiplier tubes (where the signal is transduced 
from energy to current). Next, each individual electrometer converts the current to a 
voltage signal, which is amplified and processed by the analog (A D) converter. The 
output is printed in concentration (parts per million, ppm) for all measurements. 

The operation of the PWMA is as follows. The operator first introduces an oil 
sample into the graphite tube of the Graphite Furnace Assembly. The PWMA start 
button, located on the control panel, is then depressed, and the entire analysis program 
is initiated and controlled by an "on board" microcomputer. As the graphite furnace 
executes the atomization cycles, the emission line spectra originating from the HCLs 
are focused through the graphite tube to the polychrometer where they are separated 
into individual spectral lines and are further conveyed to the exit slits. The spectral 
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lines pass through their respective exit slits where their spectral intensities are detected 
and measured by individual Photomultiplier Tubes (PMTs). Wear metals present in 
the oil sample will absorb the light associated with each metal's spectral line, thus 
reducing the intensities reaching the PMTs. The amount of light absorbed is related to 
the concentration of metal in the oil sample. 

The PWMA has several function modes. In the calibration mode, a prepared 
sample with known concentrations is used to get the absorbance data. Analysis mode is 
selected to analyze the unknown samples. The prototype PWMA prints out two types 
of data in the analysis mode, concentration data (computed by the on-board 
microprocessor) and absorbance data (in parentheses). These data can be used to 
compare analysis algorithms. 

C. TEST PROGRAM AND DATA COLLECTED 

The four prototypes were tested at their respective field evaluation sites. The test 
followed the PWMA field test burn sequence, which is in Appendix A. The graphite 
tubes, used in the graphite furnace, deteriorate gradually with use, and should be 
replaced after even - 160 sample burns. If a new graphite tube is being used, 
preconditioning of the tube is required. Laboratories were to follow the test sequence 
which included calibration samples, verification samples, correlation samples, random 
samples, reslope, and so on. Three calibration samples were prepared in MIN-L 7S0S 
oil containing 20%, 50% and 100% of full scale for all nine metals. Perkin-Elmer 
provided these samples. Three verification samples were prepared in MIN-L 7S08 oil 
containing 10%, 40% and 70% of full scale for all nine metals. Perkin-Elmer provided 
these samples, too. Correlation Samples used turbine engine oil samples. Each sample 
contained a minimum of 7 wear metals. Random samples were routine daily samples 
obtained at each evaluation site. This analysis employs only the calibration and 
verification samples. 

The concentration data printed out by the PWMA are calculated by the 
microprocessor with the use of a calibration technique. Samples with known 
concentrations are analyzed first and the absorbance data obtained are used to 
construct calibration curves. Absorbance data obtained from unknown samples are 
then compared to the calibration curves so that their concentrations can be calculated. 
Later verification samples obtained from known samples give an indication of how 
valid the calibration curve is. 
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Table 1 shows data collected from one of the field test sites. P523 means 
"Pensacola May 23rd, 1985", the place and the date the data were collected, "401" is 
the tube sequential number. CAL IA and CAL IB stands for the calibration sample 
20%. Cal II and Cal III are 50% and 100% of full scales respectively. In order to 
provide a more accurate result, a two-point average calibration was used for the field 
test, i.e., two CAL I samples were run at CAL I position; two CAL II samples at CAL 
II position; and two CAL III samples at CAL III position. For each calibration run, 
the numerical values are the absorbance data printed out (rounded to integers). The 
PWMA averaged the two samples at the same level, and stored them as the calibration 
points, but the average value was not printed out. 

TABLE 1 

ABSORBANCE DATA OF TUBE 401 



P523 401 


Element 






Fe 


Ag 


A1 


Cr 


Cu 


Mg 


Ni 


Si 


Ti 


CAL IA 




58 


53 


10 


38 


23 


44 


40 


10 


10 


CAL IB 




56 


52 


10 


34 


22 


46 


38 


8 


9 


CAL I I A 




115 


98 


19 


83 


47 


60 


92 


28 


30 


CAL I IB 




111 


100 


21 


84 


43 


61 


96 


18 


29 


CAL I I I A 




143 


129 


26 


104 


70 


72 


119 


33 


42 


CAL 1 1 IB 




145 


131 


24 


104 


71 


74 


120 


35 


44 


VERIF 10 


CONC 


9 


1 


2 


1 


4 


3 


3 


2 


5 


VERIF 10 


ABS 


28 


28 


5 


17 


12 


32 


19 


4 


11 


VERIF 40 


CONC 


44 


4 


10 


4 


18 


9 


12 


9 


7 


VERIF 40 


ABS 


105 


90 


17 


70 


42 


56 


79 


20 


20 


VERIF 70 


CONC 


70 


7 


9 


6 


28 


17 


20 


12 


10 


VERIF 70 


ABS 


129 


114 


16 


88 


56 


66 


105 


26 


29 



Part of the data for tube 401 is shown in Table 1. The complete data for tube 
401 is attached in Appendix B. VERIF 10, VERIF 40 and VERIF 70 are the 
verification samples. The data in the first row of VERIF 10 is the computed 
concentration and the second row is absorbance data. 
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Table 2 shows all the data available for this study, Langley by location tested 6 
tubes, Myrtle Beach tested 5 tubes, Pensacola tested 1 1 tubes, and Elmendorf tested 6 
tubes. All 2S tubes were used during the field test program. 

TABLE 2 

DATA OF EACH LABORATORY 



Lab. 


Tubes 


Total 


Langley 


PPM201 

PPM204 


PPM202 

PPM206 


PPM203 

PPM207 


6 


Myrtle 

Beach 


PPM302 

PPM305 


PPM303 

PPM306 


PPM304 


5 


Pensacola 


PPM401 

PPM405 

PPM409 

PPM415 


PPM403 

PPM406 

PPM412 

PPM416 


PPM404 

PPM408 

PPM413 


11 


Elmendorf 


PPM501 

PPM504 


PPM502 

PPM505 


PPM503 

PPM507 


6 



In the PWMA, three calibration standards of 20%, 50% and 100% are used to 
construct the calibration curve. The absorbance data in CAL I A through CAL II IB 
are used to compute the calibration curve coefficients. A pair of two point rational 
equations were used to fit the calibration curves, one for each of the nine elements. 
The function used is given in equation 2.1. 

C = KjA/(K 2 A-1) (2.1) 

where K [t K, are unknown coefficients, C is concentration in ppm, and A is 
absorbance data. 

The first part of the curve uses absorbance values obtained from 20% and 50% 
samples and the second part of the curve uses values obtained from 50% and 100%. 
We will explain how the concentration value of 9, the first line of the VERIF 10 in 
Table 1, was computed in the PWMA. Figure 2.1 shows the calibration technique 
used by the PWMA. In Figure 2.1, the first part of the curve for the FE is given by 
the solid line, the second part of the curve is given by the dashed line, and symbol + 
identifies the three calibration points. The verification absorption values of 129, 105 
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Figure 2.1 Calibration Technique used by the PWMA. 

and 2S are applied to the curve, and the calibrated data computed by the PWMA are 
70. 44 and 9 respectively. However, the correct values are 70, 40 and 10. Verification 
sample VFRIF 10 which has known concentration 10 produces absorbance data 28, as 
shown in Table 1. II' the absorbance data 28 was fitted to the curve from the 
absorbance axis, then its concentration should be close to 10. However, it shows some 
deviation because the curve point projected to the concentration axis makes 
concentration 9, as shown in Figure 2.1. When the observed absorbance is within the 
range of 0 to the absorbance value at 50% concentration, the first curve is applied to 
compute the concentration. For the absorbance data greater than this, the second 
curve is applied to compute the concentration. The computer selects the proper curve 
to calculate data, dependent upon the calibrating concentrations of all nine elements. 
[Ref. 2: p. 67] 

The PWMA has a multi-element capability, the 9 elements already listed in Table 
1. Iron, Silver, Aluminum. Chromium. Copper, Magnesium, Nickel, Silicon and 
Titanium. 
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The full scale capabilities for these elements are shown in Table 3. For example, 
the full scale value for Fe is 100. and is 10 for Ag. The columns labeled 20%, 50% and 
100% are CAL I. CAL II and CAL III samples respectively. The 10%. 40% and 70% 
columns are VERIF 10, VERIF 40 and VERIF 70 respectively. The concentration 
varies with the sample and the element. 

TABLE 3 

TABLE OF OIL SAMPLE MIXTURES 



METAL 


concentration (ppm) 




10% 


20% 


40% 


50% 


70% 


100% 


FE 


10 


20 


40 


50 


70 


100 


AG 


1 


2 


4 


5 


7 


10 


AL 


2. 5 


5 


10 


12. 5 


17. 5 


25 


CR 


1 


2 


4 


5 


7 


10 


CU 


4 


8 


16 


20 


28 


40 


MG 


2. 5 


5 


10 


12. 5 


LO 
i — i 


25 


NI 


3 


6 


12 


15 


21 


30 


SI 


2 


4 


8 


10 


14 


20 


TI 


2 


4 


8 


10 


14 


20 




CAL I 




CAL II 




CAL III 
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III. ALGORITHMS 



A. GENERAL 

If one increment of concentration C causes an equal increment of absorbance A 
at all levels of C, then we say that there is no curvature in the analytical curve 
expressing the relation between C and A. In reality, in the presence of stray light the 
analytical curve bends towards the concentration axis. [Ref. 3: p. 6S1] 

Calibration curvature can be controlled to an extent by the instrument design. 
The degree of curvature for iron at the 248.3-nm wavelength, for example, is partially 
controlled by the spectral bandwidth since most of the curvature can be attributed to 
overlap of nearby nonabsorbed lines in the iron spectrum. This curvature effect is 
analogous to stray light interferences observed in molecular spectroscopy. However, at 
practical instrument settings, it is not possible to eliminate this curvature completely. 

B. E. Limbeck investigated transformations of the data from AA (Atomic 
Absorption) instruments to forms that would produce more easily handled curves. 
Graphic trials on several possible transformations revealed that a plot of A C vs A or 
C A vs C produced pseudocalibrations that were much less curved than the A vs C 
calibrations, even for the extreme cases. He concluded that over a limited 
concentration range, two or three standards were sufficient to produce excellent fits for 
AA calibrations with the rational algorithm. [Ref. 4: p. 96] 

As already mentioned, the PWMA uses a rational function which uses two 
standards; to cover the full range this equation is used twice. This same function will 
be used twice in the following analysis. The PWMA on board computer uses single 
precision arithmetic and establishes its calibration curve internally. The data available 
for this study was printed out in the PWMA paper tapes, which rounded the observed 
absorbance numbers. Thus, if one uses the PWMA equation to establish a curve, and 
then reads the rounded absorbance figures back through this curve, one does not 
exactly arrive at the concentration printed out. Thus the first algorithm should give 
essentially the same results as those labeled PWMA, which were the rounded figures 
from the internal curve. The other two rational functions have three unknown 
coefficients so that they are used one time with three calibration standards. We will 
examine three rational functions and two quadratic functions, which are defined below. 

Algorithm 1 : C = A / ( Kj 4- K, A) 
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Algorithm 2 : C = ( Kj A + K 3 A 2 ) / ( K 2 A - 1) 

Algorithm 3 : C = A / ( Kj + K, A + K 3 A“ ) 

Algorithm 4 : C = Kj + K, A + A 2 

Algorithm 5 : A = Kj + K-, C + K 3 C 2 



B. ALGORITHM 1 

This is useful when a simple linear approximation gives a good fit to the plot of 
A/C versus A. The following rational function was chosen as algorithm 1: 

A / C = K, + K 2 A (3.1) 

It is clear by a simple rearrangement of the formula that the concentration is a 
function of absorbance: 



C = A / ( Kj + K 2 A ) 



(3.2) 



For brevity, the calibration algorithm using this linear rational function is 
referred to as the "rational method" and is the one actually employed in the PWMA. 
This algorithm uses a pair of two-point rational equations 3.2 to fit data for all nine 
elements. The first part of the curve uses absorbance data obtained from 20% and 
50% samples, and the second part of the curve uses data obtained from 50% and 
100% samples. 

C. ALGORITHM 2 

William B. Barnett found that the equation, 



C = 



K, A + K, A 2 

K n ! 3 - 

0 K, A - 1 



(3.3) 



was applicable in fitting the analytically useful concentration range of atomic 
absorption calibration curves. [Ref. 5: p. 829] 
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In this equation. C is the concentration. A is the observed absorbance, and the 
K/s are coefficients which must be determined during calibration. The coefficient K n is 
the reslope constant which is set equal to 1.0 during the initial standardization, which 
gives 



K, A + K 3 A 2 
K, A - 1 



(3.4) 



In equation 3.4, concentration has been expressed as a function of absorbance. This is 
the reverse of the usual representation where absorbance is considered to be a function 
of the standard concentrations. However, expressing concentration as the dependent 
variable facilitates its computation without the need for extracting roots of equations 
as others had done. Equation 3.4 can actually be used in three forms: as shown in the 
PWMA. with K, being set equal to zero, 



C = Kj A / ( K 2 A -1 ) 



(3.5) 



or with K-, and K^ both being set equal to zero. This latter case is, of course, the 
equation for a straight line which passes through the origin. 

When the number of standards matches the number of coefficients to be 
determined, it is a simple matter to determine the coefficients by the solution of 
simultaneous equations. We will employ simultaneous equations, to show how to get 
the coefficients in APL. 

Suppose 

C ( C|, c 2 . Cj ), A ( 3j, a*,, ). 

If we express both sides of equation 3.4 in terms of components and then equate 
corresponding components, we obtain the system 



K, a. + K 3 a, 2 
K 2 a ; - 1 



i = 1, 2, 3 



(3.6) 



which leads to 

K, 3j c. - c ; = Kj a ; + K 3 a. 2 



IS 



i = 1, 2, 3 



(3.7) 



i = 1, 2. 3 



(3.8) 



i.e., 

c = - Kj a ; + K-, a ; c ; - K 3 a ; 2 
Define 

K = ( K r K 2 , K 3 ). and P = || p.. || 

with 



P-, = - a., 

P i2 = a. c, 

Pi3 = - a i 2 ’ 

Then 

C = P K (3.9) 

The solution for K is 
K + C&P 

where 1 is the matrix inverse operator in APL [Ref. 6: pp. 82, 83). 

While this curve may well faithfully represent the Absorbance - Concentration 
relationship over the useful range, there are certain possible observed values where it 
does not. For example algorithm 2 will not perform well when the observed pairs (c. 
a ; ) lead to K 2 = 1/a, where aj < a < a 3 , since the denominator of equation 3.6 

becomes zero. For example, with C = ( 20, 50, 100 ), A =( 43, 117, 152 ), the 

calibration coefficients K are ( -0.4693, 0.00694, 0.003322). The fitted curve defined by 
these coefficients is shown in Figure 3.1. The horizontal line is an asymptote line at A 
equal to 144.16, which is the reciprocal of K v If the observed absorbance is slightly 
less than 144.16, the computed concentration will be a large negative number, and 
when the absorbance data is slightly greater than 144.16, then the computed 

concentration will be a large positive number, a behavior which would not be 

acceptable. 
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Figure 3.1 An Asymptote Line in Algorithm 2. 

D. ALGORITHM 3 

this algorithm was suggested by Limbeck [Ref. 4: p. 96], 

C T (3.10) 

K, + k\ A + K 3 A 2 

This is a different extension of algorithm 1, namely, A/C = + K 2 A. In 

equation 3.l<>, the quadratic term is included. To get the coefficients in APL, we again 
employ the notation of simultaneous equations. 

K, + K 2 a. + K 3 a j 2 = a./c. i = 1, 2. 3 (3.11) 

Now, let 



C = ( Cp a, c„ a 3 c 3 ) with 
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Pi = a *• 

and again C = P K, so the solution is 

K + CRP 

This algorithm can also produce a curve with an interior asymptote, like 
algorithm 2, when the denominator is zero. 

Another type of undesirable shape can occur with this case (as well as algorithm 
2). This is illustrated by taking 

C = ( 20, 50, 100 ), 

A = ( 43, 75, 152 ). 

Then, we find 

K = ( 3.632, -0.0425, 0.0001SS ). 

The resulting curve is pictured in Figure 3.2. Increased absorbance will cause lower 




Figure 3.2 Possible Observed Curve from Algorithm 3. 
estimated concentration for A > 140 with this example. 
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E. ALGORITHM 4 



A simple quadratic equation 

C = K t + K 2 A + K 3 A 2 (3.12) 

might be a reasonable alternative for a calibration curve. 

Equation 3.12 has been expressed as a function of absorbance. This is the reverse 
of the usual representation where absorbance is considered to be a function of the 
standard concentrations. However, expressing concentration as the dependent variable 
facilitates its calculation without the need for extracting roots of equations as is needed 
for algorithm 5. 

This algorithm may give a satisfactory shape, depending on the A values which 
occur during calibration burns. 

It has an undesirable shape for 

C = (20, 50. 100), 

A = (40, 130, 170), 

K = (43.3 -0.865, 0.007). 

Pictured in Figure 3.3, here decreasing A, for A < 60, gives increasing C. 

When 



C = (20, 50, 100), 

A = (40, 60, 170), 

K = (-59.3, 2.3, -0.008), 

algorithm 4 gives the curve shown in Figure 3.4, in which increasing A, for A > 140, 
leads to decreasing C. 

Note that neither of these two curves pass through the origin, since the 
polynomial equation will not necessarily give zero concentration with zero absorbance. 
In APL the coefficients K are given by 
K + C 1/1 », *0 1 2 




figure 3.3 The Curve Bending Downward at Lower Concentration by Algorithm 4. 

where C = (c ( , c v c 3 ) are the concentrations used and A = ( a^ a 2 , a 3 ) are the 
observed absorbance numbers. 

F. ALGORITHM 5 

The final curve to be considered is 

A = K, +K, C + K 3 C 2 (3.13) 

In equation 3.13, absorbance has been expressed as a function of concentration 
which is the natural representation. However, expressing absorbance as the dependent 
variable requires the solution of a quadratic equation to find concentration. Equation 
3.13 can be expressed like this, 

K 3 C 2 + K, C + K, - A = 0 (3.14) 

Solving the resulting equation for C, we obtain the formula, 
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Figure 3.4 The Curve Bending Upward at High Concentration by Algorithm 4. 



C = 



- K,±VK 2 2 - 4 K 3 ( K, - A ) 



2 K 



3 



(3.15) 



Whenever K-, 2 - 4 ( Kj - A ) is negative, there is no concentration value C for the 

given absorbance A. 

As with algorithm 4, two different undesirable shapes may occur; if the initial 
calibration data is 

C = (20, 50, 100), 

A = (40, 140, 170), 

then 



K = (-60.8, 5.7, -0.03), 

giving the curve pictured in Figure 3.5. 
If the initial data is 
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Figure 3.5 The Curve Bending Downward at High Concentration by Algorithm 5. 

C = (20, 50,100), 

A = (40, 50, 170), 

then 

K = (59.2, -1.5, 0.03), 

giving the curve in Figure 3.6. With A and C as defined earlier, the coefficients K are 
given by 

X + . *012. 



25 




Figure 3.6 The Curve Bending Upward at Lower Concentration by Algorithm 5. 



IV. ALGORITHM COMPARISON 



A. GENERAL 

We have mentioned how algorithms are used to fit the curve and each algorithm 
has delects for specific A values. Recall that the field test program included the 
concentrations: 10%, 20%, 40%, 50% 70% and 100% of full scale for each element. 
Calibrating the tube consisted of two analyses each using the 20%, 50% and 100% 
samples: these samples were used to fix the absorbance - concentration curve for all 
analyses performed until the next calibration. During the planned sequence of burns, 
each of the 10%, 40%, and 70% verification samples were analyzed twice, once just 
after the calibration curve was established and then a second time about 30 analyses 
later. Ideally, the established calibration curve should give 10%, 40% and 70% of full 
scale respectively for these samples; the different possible algorithms can be compared 
by using these 10%, 40% and 70% samples, in terms of how close the resulting 
calibration curve result is to the known concentration, for the various tubes used. 

Specifically, each algorithm can be used to establish a calibration curve using the 
observed absorbance numbers for the 20%, 50% and 100% samples. Since these curves 
are not, in general, exactly the same for all algorithms, the difference between the 
returned concentration for a 10% sample, e J0 say, can be compared with the expected 
result c lQ . We will introduce the sum of squares method to compare the algorithms. 
Thus, for a given algorithm, we will have two e )0 , two e 40 , and two e, Q values from 
which we will compute the sum of the 6 squared differences. The algorithm which has 
the smallest mean sum of squares will be selected as the rank 1, with ranks 2, 3, ..., 5 
assigned the larger mean sum of squares. The number of rank 1 and rank 2 for each 
element throughout all tubes will be considered as the comparing tool for that rank. 
The algorithm which produces the most rank 1 and rank 2 scores will be designated the 
best one. 

B. SUM OF SQUARES METHOD 

We will explain the sum of squares method by an example. For tube 401, for 
element Fe, the fitted curve is shown in Figure 4.1 where the coefficients are computed 
by algorithm 2. 

The residuals E are defined by 
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Figure 4.1 The Horizontal Deviation of Residual by Algorithm 2. 

n = c - c 

where C is a known concentration and C' is the estimated concentration produced by 
the algorithm. 

For tube 401, for Fe, the absorbance data of the verification sample A are, 

A = ( 28, 105, 129 ) 
and the known concentrations are, 

C = ( 20. 50, 100 ) 

Then we find K by algorithm 2, 

K = ( -0.3211, 0.00602, 0.00159 ) 

We can compute the estimated concentration C' from equation 3.4 with (28, 105, 129) 
in place of A. That is 

C = ( 9.315, 44.004, 67.040 ). 
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Suppose E = ( ej. e-,, ). The residual E is 

E = C - C' = ( 0.685, - 4.004. 2.96 ) 

and the sum of squares SS is defined as the sum of the squares of these residuals, 
i. e., 

3 

SS = V e 2 = 0.6S5 2 + 4.004 2 + 2.96 2 = 25.263. (4.1) 

i= 1 



Then the mean sum of squares MSS is defined as. 



MSS 



S S 

~T 



(4.2) 



giving 

MSS = 25.263 / 3 = 8.42. 

Residuals at each point are shown in Figure 4. 1 where the residual is indicated by 
the horizontal line between symbol ' x ' and 'o'. Estimated concentration C' is 
indicated by the symbol ' x ' and concentration is indicated by the symbol 'o'. 

When we are using algorithm 4 with the same data, then 

K = ( 79.237, - 1.775, 0.013 ) 

This curve is shown in Figure 4.2. Horizontal deviation is the distance between symbol 
' x ' and ’o', which is shown as a horizontal line. The estimated concentrations are 

C = ( 39.75, 36.57, 67.184 ). 

As we see, the estimated concentration at absorbance 28 is larger than at absorbance 
105, because the curve is bent downward at the lower concentration with this 
algorithm. Therefore, the residuals are larger. Also, algorithm 2 is considerably better 
than algorithm 4 for this set of data. 

Accordingly, 
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Figure 4.2 The Horizontal Deviation of Residual by Algorithm 4. 
E = ( -29.75, 3.43, 2.816 ), 

SS = (-29.75)" + 3.43" + 2.S16 2 - 735.257, 



and 



MSS = 735.257 3 = 245.086. 

For these data with these two algorithms, 

MSS by algorithm 2 = 8.42, 

MSS by algorithm 4 = 245.086, 

and algorithm 2 is better than algorithm 4. 

Estimated concentrations C' by 6 algorithms for Fe are summarized in Table 4. 
Known concentrations C, shown in the first column, were used 4 times 10%, 40% and 
70% of full scale. After the calibration curve was established, each of the 10%, 40%, 
and 70% verification sample was analyzed twice. For tube 401 this calibration was 
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done twice, so that 12 verification data are gathered in Table 4. Notice the differences 
between the estimated concentrations in column 2 and 3 in Table 4. The estimated 
concentrations by PWMA and by algorithm 1 are very close. Since algorithm PWMA 
and algorithm 1 are exactly the same as noted earlier: the PWMA concentration values 
are those printed out by the microprocessor, while the algorithm values are computed 
from the curve established from the rounded absorbance values. Since algorithms 2 
through 5 must be employed with the rounded absorbance values, rather than the 
single precision values available to the microprocessor the rank comparisons will be 
restricted to only algorithms 1 through 5. The mean square residuals will be computed 
for the PWMA, but they will not be used in comparison with the others; cases in which 
the PWMA actually gives the smallest mean square will be denoted by an asterisk(*). 
We might expect that algorithm PWMA and algorithm 1 should perform very’ 
similarly. In the last column. .00 means we can not obtain acceptable estimated 
concentrations with this algorithm. If K, 2 - 4 K 3 ( Kj- A ) was negative, the estimated 
concentration was forced to be .00 for all verification samples with this calibration. 
The estimated concentration in the first three rows of algorithm 2 and 4 are the same 
as our previous example. 

From the data in Table 4, we calculate the sum of squares and the mean sum of 
squares, which are displayed in Table 5. It says that the sum of squares by PWMA is 
100, and the mean sum of squares is 8.33, so the rank of the mean sum of squares is 
labeled with * as it is the smallest. If we don't count the PWMA in comparison, 
algorithm 3 is the best. It has the smallest mean sum of squares, 9.09: thus algorithm 
3 suggested by Limbeck is the best for the tube PPM401 for element Fe. The number 
of acceptable data points for each algorithm is given in the row labeled X; if the test 
plan was followed exactly, and the algorithm performed satisfactorily for the 
calibration samples, N would equal 12 for each case. If the algorithm performed totally 
unsatisfactorily, its N value is 0 and its mean sum of squares was set equal to 888 (see 
algorithm 5 in Table 5). 

The sum of squares, the mean sum of squares, and the rank of the algorithm for 
tube 401 for 9 elements are summarized in Table 6. The equivalent data for the other 
tubes are attached in Appendix D. 

Table 6 shows the best algorithm for each element. Algorithm 3 is the best 
algorithm for element Fe. Algorithm 1 has the smallest mean sum of squares for Ag. 
and so on. As can be seen, the best algorithm differs according to the element. For 
tube 401, no single algorithm is the best for all the elements. 
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TABLE 4 

ESTIMATED CONCENTRATION BY 6 ALGORITHMS 



Standard 


Algorithm 




PWMA 


1 


2 


3 


4 


5 


10. 00 


9. 00 


8. 87 


9. 31 


13. 44 


39. 75 


. 00 


40. 00 


44. 00 


44. 79 


44. 00 


45. 66 


36. 57 


. 00 


70. 00 


70. 00 


70. 23 


67. 04 


67. 96 


67. 18 


. 00 


10. 00 


13. 00 


13. 21 


13. 57 


17. 64 


29. 09 


. 00 


40. 00 


42. 00 


42. 32 


41. 36 


43. 08 


32. 93 


. 00 


70. 00 


70. 00 


70. 23 


67. 04 


67. 96 


67. 18 


. 00 


10. 00 


8. 00 


8. 00 


8. 34 


17. 04 


53. 22 


. 00 


40. 00 


43. 00 


43. 05 


42. 64 


60. 82 


36. 86 


. 00 


70. 00 


74. 00 


75. 32 


74. 43 


88. 77 


83. 90 


. 00 


10. 00 


6. 00 


6. 07 


6. 39 


13. 58 


66. 29 


. 00 


40. 00 


36. 00 


35. 69 


35. 19 


52. 89 


26. 16 


. 00 


70. 00 


73. 00 


73. 52 


72. 64 


87. 41 


81. 58 


. 00 



TABLE 5 

ALGORITHM COMPARISON FOR FE 





Algorithm 




PWMA 


1 


2 


3 


4 


5 


SS( Fe) 


100. 00 


128. 08 


121. 12 


109. 09 


3119. 08 


. 00 


N 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


. 00 


MSS 


8. 33 


10. 67 


10. 09 


9. 09 


259. 92 


888. 00 


RANK 


* 


3 


2 


1 


4 


5 



In Chapter II we mentioned that Langley used 6 tubes, Myrtle Beach used 5 
tubes, Pensacola used 1 1 tubes, and Elmendorf used 6 tubes. Due to missing data, we 
excluded several tubes in our analyses. For example, tube PPM406 doesn't have 
verification sample 70 data. For similar reasons we discard 1 tube from Langley, 1 



TABLE 6 

SUMMARY OF THE RANKS FOR TUBE 401 





Algorithm 




PWMA 


1 


2 


3 


4 


5 


SS( Fe) 


100. 00 


128. 03 


121. 08 


109. 09 


3119. 08 


. 00 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


. 00 


MSS 


8. 33 


10. 67 


10. 09 


9. 09 


259. 92 


888. 00 


RANK 


* 


3 


2 


1 


4 


5 


SS( Ag) 


2. 00 


. 85 


1. 70 


1. 48 


12. 34 


. 67 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


6. 00 


MSS 


. 17 


. 07 


. 14 


. 12 


1. 03 


. 11 


RANK 




1 


4 


3 


5 


2 


SS( Al) 


268. 00 


317. 88 


300. 21 


297. 57 


429. 29 


. 00 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


. 00 


MSS 


22. 33 


26. 49 


25. 02 


24. 80 


35. 77 


888. 00 


RANK 


* 


3 


2 


1 


4 


5 


SS( Cr) 


16. 00 


17. 29 


17. 16 


17. 10 


39. 32 


. 00 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


. 00 


MSS 


1. 33 


1. 44 


1. 43 


1. 42 


3. 28 


888. 00 


RANK 


* 


3 


2 


1 


4 


5 


SS( Cu) 


106. 00 


102. 05 


100. 06 


99. 74 


100. 09 


107. 97 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


MSS 


8. 83 


8. 50 


8. 34 


8. 31 


8. 34 


9. 00 


RANK 




4 


2 


1 


3 


5 


SS( Mg) 


91. 00 


170. 50 


156. 21 


148. 12 


179. 10 


. 00 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


. 00 


MSS 


7. 58 


14. 21 


13. 02 


12. 34 


14. 92 


888. 00 


RANK 


* 


3 


2 


1 


4 


5 


SS(Ni) 


38. 00 


34. 87 


26. 02 


31. 70 


192. 41 


. 00 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


. 00 


MSS 


3. 17 


2. 91 


2. 17 


2. 64 


16. 03 


888. 00 


RANK 




3 


1 


2 


4 


5 


SS( Si ) 


162. 00 


149. 09 


141. 29 


166. 29 


171. 89 


103. 04 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


6. 00 


MSS 


13. 50 


12. 42 


11. 77 


13. 86 


14. 32 


17. 17 


RANK 




2 


1 


3 


4 


5 


SS( Ti ) 


135. 00 


142. 77 


364. 12 


141. 42 


134. 41 


141. 47 


n 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


12. 00 


MSS 


11. 25 


11. 90 


30. 34 


11. 79 


11. 20 


11. 79 


RANK 




4 


5 


2 


1 


3 



tube from Myrtle Beach, 2 tubes from Pensacola, leaving 24 tubes in all. The summary' 
of 24 tubes for each element is shown in Table 7 through Table 15. We may say that 
algorithms which have the most rank 1 or rank 2 scores but the least rank 4 or rank 5 
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TABLE 7 

FLANKS OF 24 TUBES FOR FE 



Algo. 


rank 




* 


1 


2 


3 


4 


5 


PWMA 


10 












1 




7 


1 


11 


5 


0 


2 




7 


8 


3 


0 


5 


3 




5 


12 


6 


1 


0 


4 




1 


0 


0 


12 


11 


5 




4 


3 


3 


6 


8 



TABLE S 

RANKS OF 24 TUBES FOR AG 



Algo. 


rank 




* 


1 


2 


3 


4 


5 


PWMA 


6 












1 




7 


3 


4 


8 


2 


2 




7 


8 


5 


2 


2 


3 




0 


7 


11 


6 


0 


4 




1 


3 


1 


5 


13 


5 




9 


3 


3 


3 


6 



scores are good ones, and algorithms which have the least rank 1 or rank 2 scores but 
the most rank 4 or rank 5 scores are not good ones. The entries in these tables are the 
counts of the ranks. For example, in Table 7 algorithm 1 for Fe has 7 rank I s, 1 rank 
2, 11 rank 3's, 5 rank 4's, and 0 rank 5 out of 24. For Fe, the PWMA internal 
algorithm, with rounded concentrations, had the smallest mean sum of squares 10 time 
out of 24. In Table 7, we conclude that algorithm 1 is better than algorithms 2, 3, 4. 
and 5. Algorithm 2 has 7 rank l's but also 5 rank 5's. Therefore, algorithm 2 is not as 
good as algorithm 1. If algorithm 2 did not have those undesirable cases with 
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TABLE 9 

RANKS OF 24 TUBES FOR AL 



Algo. 


rank 




X 


1 


2 


3 


4 


5 


PWMA 


12 












1 




2 


8 


6 


4 


4 


2 




10 


4 


7 


3 


0 


3 




6 


8 


9 


1 


0 


4 




3 


1 


1 


10 


9 


5 




3 


3 


1 


6 


11 



TABLE 10 

RANKS OF 24 TUBES FOR CR 



Algo. 


rank 




* 


1 


2 


3 


4 


5 


PWMA 


9 












1 




9 


2 


7 


5 


1 


2 




4 


6 


10 


2 


2 


3 




3 


15 


4 


2 


0 


4 




0 


0 


2 


11 


11 


5 




8 


1 


0 


4 


10 



denominator zero, then it would be the best. Also algorithm 4 has 1 1 rank 5's because 
it has a lot of deviations for two extreme cases as we saw in the previous Chapter. 
Algorithm 5 is almost the same as algorithm 4. Algorithm 3, with 5 rank l's and 12 
rank 2's. is acceptable but not as good as algorithm 1. We conclude that algorithm 1 
is the best for element Fe. 

For element Ag, even though algorithm 5 has 9 rank l's, it also has 6 rank 5's. 
Therefore, we can not select algorithm 5 as the best one. Algorithm 2, or 1 are 
acceptable for element Ag. 
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TABLE 11 

RANKS OF 24 TUBES FOR CU 



Algo. 


rank 




* 


1 


2 


3 


4 


5 


PWMA 


15 












1 




12 


1 


1 


9 


1 


2 




4 


8 


5 


6 


1 


3 




1 


7 


9 


0 


2 


4 




2 


2 


3 


4 


13 


5 




5 


6 


1 


5 


7 



TABLE 12 

RANKS OF 24 TUBES FOR MG 



Algo. 


rank 




* 


1 


2 


3 


4 


5 


PWMA 


12 












1 




12 


5 


6 


0 


1 


2 




4 


7 


11 


2 


0 


3 




6 


9 


4 


5 


0 


4 




0 


2 


3 


14 


5 


5 




2 


1 


0 


3 


18 



For Al, algorithm 2 dominates the other algorithms, since it has 10 rank l's. 
Algorithm 3 is the second algorithm. Algorithm 3 could be useable even though it has 
fewer rank 1 scores. 

For Cr, algorithm 4 is not acceptable, for it has 22 tubes rank 4's or 5's. 
Algorithm 1 is the best for element Cr. 

For Cu, Mg, and Ni, algorithm 1 dominates the others. For Si, algorithm 2 or 3 
are the best. For Ti, algorithm 2 is good. Algorithm 1 and 4 are acceptable too. 
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TABLE 13 

RANKS OF 24 TUBES FOR NT 



Algo. 


rank 




* 


1 


2 


3 


4 


5 


PWMA 


7 












1 




13 


2 


8 


1 


0 


2 




4 


2 


2 


5 


11 


3 




4 


18 


2 


0 


0 


4 




1 


1 


3 


16 


3 


5 




2 


1 


9 


2 


10 



TABLE 14 

RANKS OF 24 TUBES FOR SI 



Algo. 


rank 




* 


1 


2 


3 


4 


5 


PWMA 


11 












1 




3 


7 


8 


3 


2 


2 




8 


7 


2 


3 


4 


3 




7 


2 


3 


11 


1 


4 




1 


7 


9 


4 


3 


5 




4 


1 


2 


3 


13 



Table 16 summarizes the results of the previous tables by listing the best and 
alternate algorithm for each element. We can say that algorithm I provides the best fit 
for elements Fe, Cr, Cu, Mg, and Ni. and algorithm 2 provides the best fit for elements 
Ag, Al, Si. and Ti. 
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TABLE 15 

RANKS OF 24 TUBES FOR TI 



Algo. 


rank 




k 


1 


2 


3 


4 


5 


PWMA 


11 












1 




5 


9 


3 


5 


2 


2 




7 


3 


2 


8 


4 


3 




4 


1 


6 


7 


6 


4 




5 


5 


11 


2 


1 


5 




3 


5 


3 


2 


11 



TABLE 16 

THE BEST ALGORITHM FOR EACH ELEMENT 





Element 




Fe 


Ag 


A1 


Cr 


Cu 


Mg 


Ni 


Si 


Ti 


First 


1 


2 


2 


1 


1 


1 


i 


2 


2 


Second 


3 


1 


3 


3 


2 


3 


3 


3 


1 / 4 



3S 



V. CONCLUSIONS AND RECOMMENDATIONS 



There is no single best algorithm for all nine elements, since the degree of 
curvature is different according to the element. The PWMA didn't consider different 
algorithms for different elements. The computation in the PWMA used only by one 
algorithm for all nine elements. If we employ different algorithms for different 
elements in an analyzer, it will be better because each element has different 
characterises. We recommend that one may use algorithm 1 for elements Fe, Cr, Cu, 
Mg. and Xi and algorithm 2 for elements Ag, Al, Si. and Ti to obtain better estimation 
in computing concentration. The microprocessor can easily check, for the existence of 
an asymptote. If one occurs, use algorithm 1. 
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APPENDIX A 

PWMA FIELD TEST BURN SEQUENCE 



Day 1 

1. Tubs Condition 

2. Auto Zero - Operate in AZ mode. If the A Z failed* repeat the AZ and 

number the run as 2A. If the AZ failed twice in a row, 
replace the tube. 

3. CAL la - Operate in CAL I mode. 

4. CAL lb - Operate in CAL I mode. 

5. CAL Ila - Operate in CAL II mode. 

6. CAL lib - Operate in CAL II mode. 

7. CAL Ilia - Operate in CAL III mode. 

8. CAL Illb - Operate in CAL III mode. If the calibration failed, 

repeat the calibration and number the calibration runs as 
3A> ...., 8A, accordingly. 

9. Blank - Run a burn with no sample j operate in Analysis mode. From 

now on, all the samples will be run in Analysis mode except 
Reslope. 

10. Verification Sample # 1: 102 

11. Verification Sample U 2: 40 2 

12. Verification Sample # 3: 702 

13 Blank 

14-18. Correlation sample # 1 - Repeat the same sample five times, this 

procedure applies for all the correlation sample. 

19-23 . Correlation sample # 2 
24-28 . Correlation sample # 3 
29. Blank 

30-37. Random samples - These may be real engine oil sample. If there are 

not enough real engine oil samples, repeat one sample 
as many times as necessary. ( 8 samples total ) 

38. 10 2 sample 

39. 402 sample 

40 70 2 sample 

Day 2 

41. Blank 

42. Reslope A - Operate in RSLP mode. If the reslope A failed, disregard 

and continue the sequence. 

43. Reslope B - Operate in RSLP mode. If the reslope B failed, disregard 

and continue the sequence. 

44. 10 2 sample 

45. 402 sample 

46. 70 2 sample 

47-51 . Correlation sample # 4 
52-56 .Correlation sample # 5 
57. Blank 

58-63. Random samples ( 6 samples total ) 

64. CAL I check sample - Operate in Analysis mode. 

65. CAL II check sample - Operate in Analysis mode. 

66. CAL III check sample - Operate in Analysis mode. 

67. Blank 

68-77. Random samples ( 10 samples total ) 



78. 


102 


sample 


79. 


402 


sample 


80. 


702 


sample 


Day 


3 




81. 


Blank 


82. 


CAL 


la. 


83. 


CAL 


Ib. 



40 



84 . CAL Ila. 

85. CAL lib. 

86. CAL Ilia. 

87. CAL Illb. 

88. Blank 

89. 10/ sample 

90. 40/ sample 

91. 70/ sample 

92-96 . Correlation sample # 6 
97-101 .Correlation sample # 7 
102. Blank 

103-117. Random samples ( 15 samples total ) 

118. 10/ sample 

119. 40/ sample 

120. 70/ sample 

Day 4 

121. Blank 

122. Reslope A 

123. Reslope B 

124. 10/ sample 

125. 40/ sample 

126. 70/ sample 

127-131 .Correlation sample # 8 
132-136 .Correia tion sample # 9 
137. Blank 

138-157 . Random samples ( 20 samples total ) 

158. CAL I check sample 

159. CAL II check sample 

160. CAL III check sample 
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APPENDIX B 
DATA PPM401 



There are several data, \v 

PPM901 



P523 901 
CAL IA 


P 


3 


58 


CAL IB 


P 


9 


56 


CAL IIA 


P 


5 


115 


CAL IIB 


P 


6 


111 


CAL IIIA 


P 


7 


193 


CAL IIIB 


P 


8 


195 


VERIF 


10 


P 


10 


9 


VERIF 


10 


ABSORBANCES 




28 


VERIF 


90 


P 


11 


99 


VERIF 


90 


ABSORBANCES 




105 


VERIF 


70 


P 


12 


70 


VERIF 


70 


ABSORBANCES 




129 


VERIF 


10 


P 


38 


13 


VERIF 


10 


ABSORBANCES 




90 


VERIF 


90 


P 


39 


92 


VERIF 


90 


ABSORBANCES 




101 


VERIF 


70 


P 


90 


70 


VERIF 


70 


ABSORBANCES 




129 


P529 901 
RESLOPE 1 


P 


92 


117 


RESLOPE 2 


P 


93 


119 


VERIF 


10 


P 


99 


12 


VERIF 


10 


ABSORBANCES 




36 


VERIF 


90 


SI >LIMIT P 


95 


96 


VERIF 


90 


ABSORBANCES 




107 


VERIF 


70 


SI >LIMIT P 


96 


79 


VERIF 


70 


ABSORBANCES 




131 


CAL I 




P 


69 


23 


CAL I 




ABSORBANCES 




65 


CAL II 


P 


65 


50 


CAL II 


ABSORBANCES 




113 


CAL III 


CR>NI,SI,P 


66 


103 


CAL III 


ABSORBANCES 




199 


VERIF 


10 


SI>LIMIT?P 


78 


9 


VERIF 


10 


ABSORBANCES 




28 


VERIF 


90 


P 


79 


91 


VERIF 


90 


ABSORBANCES 




99 


VERIF 


70 


SI >LIMIT P 


80 


73 


VERIF 


70 


ABSORBANCES 




131 



just one of them for references. 



10 


38 


23 


99 


90 


10 


10 


10 


39 


22 


96 


38 


8 


9 


19 


83 


97 


60 


92 


28 


30 


21 


89 


93 


61 


96 


18 


29 


26 


109 


70 


72 


119 


33 


92 


29 


109 


71 


79 


120 


35 


99 


2 


1 


9 


3 


3 


2 


5 


5 


17 


12 


32 


19 


9 


11 


10 


9 


18 


9 


12 


9 


7 


17 


70 


92 


56 


79 


20 


20 


9 


6 


28 


17 


20 


12 


10 


16 


88 


56 


66 


105 


26 


29 


8 


1 


8 


2 


9 


9 


8 


15 


25 


22 


30 


25 


21 


21 


7 


5 


16 


8 


12 


15 


7 


19 


77 


39 


59 


76 


30 


18 


11 


7 


26 


19 


21 


12 


12 


18 


99 


59 


68 


107 


26 


39 


12 


72 


51 


55 


80 


37 


19 


13 


69 


51 


59 


80 


39 


17 


1 


1 


9 


2 


3 


9 


3 


3 


17 


13 


28 


19 


10 


7 


8 


9 


20 


8 


11 


25 


8 


15 


79 


95 


53 


70 


92 


21 


10 


7 


31 


19 


17 


39 


10 


16 


95 


61 


63 


99 


67 


30 


3 


2 


8 


9 


5 


8 


9 


7 


38 


23 


92 


39 


17 


10 


9 


6 


20 


19 


19 


19 


9 


16 


88 


99 


61 


89 


28 


25 


28 


11 


35 


25 


31 


30 


23 


28 


110 


65 


79 


129 


50 


50 


2 


1 


9 


2 


2 


0 


3 


5 


21 


12 


30 


15 


117 


6 


11 


5 


15 


11 


12 


17 


8 


18 


79 


37 


59 


76 


31 


23 


22 


7 


25 


23 


21 


23 


13 


29 


95 


53 


72 


107 


39 


39 



e list 

53 

52 

98 

100 

129 

131 

1 

28 

9 

90 

7 

m 

i 

28 

9 

89 

7 

115 

101 

96 

1 

29 

9 

91 

9 

129 

2 

52 

5 

98 

10 

130 

1 

29 

9 

80 

7 

119 
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APPENDIX C 

PROGRAMMING LIST 



Cl] 

[ 2 ] 

[3] 

14 ' 

[5] 

[ 6 ] 
[71 
[ 8 ] 
[9] 
[10 
[11 

[12 

[13 

[14 

[15 

[16 

[17 

[18 

[19 

[20 

[21 

[22: 

[23 

[24 

[25 

[26 

[27 

[28 

[29 

[30 

[31 

[32 

[33 

[34 

[35 

[36 

[37 

[38 

[39 

[40 

[41 

[42 

[43 

[44 

[45 

[46 

[47 

[48 

[49 

[50 



i 5 ] 5*2 

01 ), 



[0 .1] (A 12 3] * 50 100 



V COEF 1 

n 8/12/86 PROGRAM TO FIT BILL ' S 2 POINT RATIONAL EQUATION . CHECK OUT 
fl VERIFICATION READINGS AND THE VALUES AFTER THE RESLOPE AS WELL . 

H+ 10. 10. 250. 10. 40. 250. 30. 20. 2 

E-6-20 + 5X i 9xJ-e-l 

Tl + 8 3 6 1 pO 

R1+ 3 0 9 6 1 pO 

Q\*- 6 0 J 6 1 pO 

'FOR TUBE 401 ' 

'MODEL A/ (K1+K2A ) 1 

] LOOP:A++J 3 2 p(<£> ,PPM401[l+x6;P[I]+t5] 

] D^(UCl 2]^ 20 50 xfl[J] )|J[l 2]o,* - ' 

1P[I] ) A[_2 3]°.* 0 1 

fl ORIGINAL COEFFICIENTS FOR THE TWO CURVES . 

D+& .PPM401 [7 + tl2;E[T]+i5] 

O+Dl 2xi6]*(, (0[2xi6]°.* 0 1) + . *§K )[1 3 6 7 9 12] 

U 1+ 83^61 p 0_ 

S-<-6 0 56 1 pD[ l + 2xt6] 

Ql+Ql, [l] Q 

^(60 56 1pD[ l + 2x i 6] ) , £/l 
Tl+Tl , [1] U1 

R+ 3 0 $ 6 1 p 10 40 70 xfl[ (E[I]*5 )-4] 

Rl*-Rl [1] R 

' for Element : » , ( el [ i ] ) , 1 .prior to re slope » 

U+U+C3 0 5 6 1 p 10 40 70 xfl[ (E[J1*5 )-4] ) ,U 
O+/U ,PPM401[20 + x2 ;E[I] + x5] )*2 
C^5 0xfl[I]f (2_pC)-S-(tf+.xC* 0 1) 

n GIVES MULTIPLICATIVE CORRECTION FACTORS FOR THE TWO LINES 
D+* .PPM401 [( 22+16) ,34+t6 ;E[J]+x 5] 

O+C L6p 1 1 2]xZ?[2x 1 6]*(,(D[2x 1 6]o.*0 l)+.x$tf)[l 36 79 12] 

1 AFTER THE RE SLOPE ' 

U+( 6 0561 pL[ 1 + 2 x x 6] ) , 8 3 5 6 1 pO 
D+b.PPM 401 [ (28 + i6 ) ;P[I]+i5] 

O+CL 1 1 2 2]x£)[2x 1 2 2 3 ] * ( , (I [ 2x 1 2 3 ] ° , 

CM 3 0 +£/),[l](((6 0541 pi [ 1 3 3 5]), 8 
OCM 3 0 5 10 1 p 10 40 70 20 50 50 100 10 40 70 xff[(E 
*( (1*1 + 1 ):£9)/L00P 
1*1 

LOOT A*-+ / 3 2 p (<l> ,PPM401[41 + x6 ;E[I] +\5] 

" ‘ ' ~ " ' @A[12' 



* 0 1 ) + . [ 1 
3 5 4 1 pO) .11] 



3 4 6] 

3 0 + CO 
5)-4] ),U 



•■15 

2] o.* O' 1), [6.1] (AZ2 3]* 50 100 



)*2 



K*( U[1 2] * 20 50 XHLI1 ) 

SP[I] ) AL2 3]°.* 0 1 
P*<i> .PPM401 [47 + 112 ;E[I] +i 5] 

0*L[2x i6] t ( , (P[2xi6]°.* 0 l)+.x §K ) [1 3 6 7 9 12] 
C/1* 8 3 5 6 1 p 0 _ 

Q* 6 0 5 6 1 pD[ 1 + 2x x 6] 

Qlt$l» [1] Q 

^(6 0 5 6 1 pD[ l + 2x x 6] ) ,U1 
Tl+Tl , [1] ill 

^ 3 0 5 6 1 p 10 40 70 xff[(£[J]45)-4] 

Rl+Rl, [1] R 

' FOR ELEMENT : ' , (EL [I; ] ) , ' , PRIOR TO RESLOPE ' 
□*[/*( 3 0®61pl04070 xfl[(E[J]*5)-4]),£/ 

*( (1*1+1 )£9 ) /LOOT 



[ 1 ] 

[ 2 ] 

[3] 

" 4 " 

: 5 : 

[ 6 ] 



V COEFR 

B 8/12/86 PP0CP4A2 PO PIP PILL ' S 2 POIiVP RATIONAL EQUATION . CHECK OUT 
B 7EPIEICAPI0A/ PP4LIA7GS d/VZ? PPE 74LC/ES 7LPPEP THE RESLOPE AS WELL . 

P* 1 0.1 0.25 0.1 0.4 0.25 0.3 0.2 0.2 
P*20+5x x 9xJ*l 
P2* 8 3 5 6 1 pO 
1 FOR TUBE 401 ' 
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'.71 

! 8 ] 

!9] 

:io 
:n 
[12 
:i3 
[14 
[15 
[16 
[17 
:i8 
:i9 
! 20 
!21 
!22 

[23 

[24 

[25 

[26 

[27 

[28 

[29 

[30 

[31 

[32 

[33 

[34 

[35 

[36 

[37 

[38 

[39 

[40 

[41 

[42 

[43 

[44 

[45 

[46 

[47 

[48 



LOOP:A*+/ 3 2 p ($ , PPM401 [l+x6;£[J]+\5] ) + 2 
X+((A«.x 1 l),A*2)t(2i) 50 100 xfl[I] )o.* l o 1 
Q*K* 3 1 p (B/f )+. x 111 

A ORIGINAL COEFFICIENTS FOR THE TWO CURVES. 

ZK<j>,PPM401[7 + xl2:£[I]+i5] 

0«-((K[ls]xZ>[2xi6j)+X[3;]xZ?[2xi6]xD[2xi6] )*((K[2;]x0[2x X 6] )-l) 
U2* 8 3 9 6 1 pO 

0-^(7 1 5 6 1 p0[ 1 + 2 x x 6 ] ) , 02 
T2+T2, [1] 02 

' FOR ELEMENT : * , ( [X ; D ) , 1 , PRIOR TO RESLOPE ' 

l>0«-(4 1 $ 6 1 p 10 40 70 xff [ (£[1] +5 )-4] ) ,U 

0+/U,PPM401[20 + i2;£[I]+x5] ) + 2 

C-6-50xp[J]*(X[l;]xC+£[3:]xCxC)+((£[2;]xC)-l) 

n GiT£S MULTIPLICATIVE CORRECTION FACTORS FOR THE TWO LINES 

D*& , PPM401 [ ( 22 + x 6 ) ,34+x6:£[J] + i 5 ] 

0**-CL6pl]x(X[l;]xZJ[2xi6]+X[3 j]xZ?[2x X 6]xU[2xi6] )+ ( Of [2 5 ] x£ [2 
6] )-l) 

' AFTER THE RESL.0PE ' 

U*( 7 1®61 gZ>[ 1 + 2 x x 6 ] ) , 8 3 5 6 1 p0 
£•«-$ .PPM401 i(28+i6);£[J]+i5] 

0^C[4pl]x(X[l;]xC[2 446]+£[3;]x£)[2446]x£[2446])*((£[ 

2:1 D[2 4 4 6] )-l ) 

•*( ( + /0S2OOO ) = 1 )/£££ 

(/♦( 3 0 +l/),[l] (((7 1 9 4 1 pfl[l 3 3 5]), 83541 p0),[l] 30+0) 

D*U* ( 3 0 5 10 1 p 10 40 70 20 50 50 100 10 40 70 x£[ (£[1] + 5 )-4] ),0 

•>( (1*1 + 1 )£9)/L00P 

1*1 

LOOT : A*+/ 3 2 p (a ,PPM401 [41 +x 6 ; £[J] +x 5] )+2 
£■<■ ( (A 0 . x 1 - l),A*2 )+(20 50 100 x£[J] )o.* 1 0 1 
3 1 n ( [|]j£ )+ x 

a ORIGINAL COEFFICIENTS FOR THE TWO CURVES . 

D*$. .PPM4 01 [47+xl2:£[J]+x5] 

0+(U[l;]xZ?[2xx6] )+£[3 ;]x£[2xx6]xD[2xx 6] )+ ( (£ [ 2 ; ] x£ [ 2 x x 6] )-l) 
U 2* 8 3 5 6 1 p 0 

U*(7 1 5 6 1 p£[ 1+2 x x 6] ) ,02 
T2*T2, [1] 02 

' £0£ ELEMENT :',(££[!;]),' .PRIOR TO RE SLOPE ' 

□*0*(4 1 5 6 1 p 10 40 70 x0[(P[J]*5)-4] ) , 0 
-»■ ( (P-«-J+l )£ 9 ) /LOOT 
*0 

ERR :_0-<-4pO 

U*( 3 0 +0 ) , [ 1 ] (((7 1 5 4 1 pD [1 3 3 5] ) , 8 3 5 4 1 p0),[l] 3 0 +0) 
□♦0«-(4 1 5 10 1 p 10 40 70 20 50 50 100 10 40 70 xff[ (£[1] + 5 )-4] ),0 
•*■( (J-«-J+l )£9 ~) /LOOP 



[1] 
[ 2 ] 
[3 
4 
[5 
:6 

'.81 

!9" 

:io] 

[ii] 

:i2] 

:i3] 

:i4] 

[15] 

[16] 
1171 
[18] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 
723] 



V COEF 3 

fl 9/06/86 PROGRAM TO FIT BILL ' S 2 POINT RATIONAL EQUATION , CHECKOUT 
fl VERIFICATION READINGS AND THE VALUES AFTER THE RESLOPE AS WELL . 
fl MODEL C - Ai (Cl+C2x^+C3x^x^ ) 

• MODEL C = Ai (Cl+C2xA+C3*AxA ) ' 

H* 1 0.1 0.25 0.1 0.4 0.25 0.3 0.2 0.2 
T 3* 8 3 5 6 1 p 0 
£•«- 2 0 + 5 x x 9 x I* 1 
' FOR TUBE 401 ' 

L00P:A*+/ 3 2 p ($ ,PPM201[l+x6 ;£[I] +x5] )+2 
□«-£<- U + O 2 0 50 100 xff[j] )BA® .* 0 1 2 
fl ORIGINAL COEFFICIENTS FOR THE TWO CURVES. 

D*0> .PPM401 [7 + xl2;£[J] + x5] 

0+D [2xx6] + (,(0[2xx6]o.* 0 1 2 )+ . x$£) 

0 3* 8 3 5 6 1 p0 

U * (60561 p0 [~l + 2x x6] ) ,03 
T3+T3 [1] u3 

1 for Element: ',(££[!;]),' .prior to re slope' 

□+0*(3 0 5 6 1 p 10 40 70 x0[ (e[J]+5 ) - 4] ),0 
C*+ / ( $ ,PPM40 1 [20+x2;£[J]+x5])+2 

0 1 2 ) 

fl CIT£S MULTIPLICATIVE CORRECTION FACTORS FOR THE TWO LINES 
D+S> .PPM4 01 [ (22 + x6 ) ,34 + x6 ;£[J] +x 5] 

0-e-CL6pl]x£)[2xx6] + (, (P[2xx6]°.*012)+.xfS)£) 
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: 24 ] 
; 2 5 D 
[26] 
[2 7] 
[28 
! 29 

; 3 o ] 

131 
[32 
; 3 3 
! 3 4 
[3 5 

; 3 6 ] 
: 3 7 ~ 

[38 
; 3 9 
[40 
[41 
[42 



:i] 

: 2 ] 
[3] 
'4] 

: 5 ] 
[6] 

[7] 

[ 8 ] 
[9] 
! 10 ] 
:n“ 
: 12 

:i3] 
:i4 _ 
:i5 
: i 6 ] 
: 1 7 ~ 
: 1 8 
:i9 

: 20 
[21 
'.22 
'.23 
; 24 
[25 
: 26 
: 2 7 
[28 
! 29 
[30 
[31 
; 3 2 
! 3 3 
[34] 
[35~ 
: 3 6 
[3 7] 
! 3 8 D 
13 9] 

[ 40 ] 

[ 41 ] 
[42 



] 3 0 40) 

5 ) -4] ),0 



' AFTER THE RE SLOPE ' 

0-5-( 6 0 561 p0[ l + 2x t 6] ) , 8 3 * 6 1 p0 
0«-4 .PPM401 [ (28+i6);F[J]+i5] 

0-5-C£4pl]x0[2x 1 2 2 3]*(,(0[2x 1 2 2 31°. *0 1 2)+.x$tf) 

0-5- ( 3 0 +£/) , [1] (( (6 0 5 4 1 pD[l 3 3 5] ), 8 3 5 4 1 pO) , [1 
□-5-0-5-(3 0 5 10 lp 10 40 70 20 50 50 100 10 40 70 xHZiEZIl* 

-*( (1*1+1 )£9 ) /LOOP 
1 + 1 

LOOT:A++/ 3 2 p ($,PPM201[41 + i6;ff[J]+i5] )t 2 
□-5-F-5-U + C* 20 50 100 xp[J] )i4o . * 0 1 2 
p ORIGINAL COEFFICIENTS FOR THE TWO CURVES . 

0*4 .PPM4 01 [47 + 1 12 ;E[J] +1 5] 

0*0[2xi6]+(, (C[2xi6]f.*01 2)+.x§F) 

03* 8 3 5 6 1 pO 

0* ( 6 0 $ 6 1 pP[~l+2xi6] ) , 03 
T3*-T3 [1] 0 3 

' for Element : « , {el z j ; ] ) , « .prior to re slope « 

□*0*(3 0 5 6 1 p 10 40 70 xtf [ (£[/] 5-5 )-4] ) ,0 
+((I+I+l)Z9)/LOOT 

V C0FF4 

p 9/09/86 PROGRAM TO FIT BILL ' 5 2 POINT RATIONAL EQUATION . CHECKOUT 
p t'FPJFICdZ’iW PF40J00S 400 PPF 7400FS 4FPFF F0F RESLOPE AS WELL . 
p MODEL C = Cl+C2x4+C3x4x4 
'MODEL C = Cl+C2xA+C3xA*A ' 

0* 10. 10. 250. 10. 40. 250. 30. 20. 2 
F*20+5x i 9x J-t-i 
P4* 8 3 $ 6 1 pO 
' FOR TUBE 401 * 

LOOP:A++/ 3 2 p ($ , PPM501 Cl + i6;FCJ]+i5] )*2 
□*F*(C* 20 50 100 xfl[I] )@4°.* 0 12 
p ORIGINAL COEFFICIENTS FOR THE TWO CURVES. 

0*4 .PPM401 [7+i 12 ;FCJ] + t 5] 

0*, (P[2xi6]°.* 0 1 2)+.*K 
U 44- 8 3 5 6 1 p0 
2’4-6-2’4 , [1] 04 

0-s- (6 0 * 6 1 pP["l+2xi6] ) ,04 

' FOR ELEMENT : « , (EL II ; ] ) , « .PRIOR TO RESLOPE ' 

□*0*(3 0 5 6 1 p 10 40 70 x0[(£[J 1+5 )-4] ),0 
C*+/(4,PPM401[20 + i2;E[.n+i5] )*2 
C*(5Ox0[J] )i(K+.*C* 0 12) 

p GIVES MULTIPLICATIVE CORRECTION FACTORS FOR THE TWO LINES 
0*4 , PPM4 01[(22 + i6),34 + i6;F CP] -K5] 

0*Cx (,(C[2xi6]o.* 0 1 2 )+. xtf) 

' ^FPFP THE RE SLOPE ' 

0* (6 0 5 6 1 qDZ 1 + 2xi 6] ) , 8 3 5 6 1 pO 
, PPM401 1 (28 + x6);F[J]+i5] 

0*Cx ( , (0[2x 1 2 2 3]°.* 0 1 2 ) + . x# ) 

0* ( - 3 0 40), [1] (((6 0 $ 4 1 pP[l 3 3 5] ), 83*41 p0),[l] 30 40) 
□*0*(3 0 * 10 1 p 10 40 70 20 50 50 100 10 40 70 xtf [ (E [I ] 4 5 ) -4] ),0 
*((J*J+1)£9)/000P 
J*1 

LOOT:A++/ 3 2 p (4 ,PPA7501 [41 + 1 6;F[J]+x5] )*2 
□*F*(0* 20 50 100 x0[J] )fflA®.* 012 
p ORIGINAL COEFFICIENTS FOR THE TWO CURVES. 

0-«-4 .PPM401 [47 + il2;F[J]+i5] 

0-5- , v0[2xi6]°.* 0 1 2) + .*K 
04-5- 8 3 * 6 1 p0 
P4-5-P4 , [1] 04 

0-5- (6 0 * 6 1 p0[ 1 + 2 x x 3 ] ) 04 
• FOR ELEMENT (EL II ; ] ) , ' .PRIOR TO RE SLOPE ' 

□-5-0-5- (3 0 * 6 1 p 10 40 7 0 Xff [ (EZI1 4 5 )-4] ) ,0 
*( (J-6-J+1 )<.9)/LOOT 



:i 
[ 2 ] 
: s ~ 

[4 



V COEF 5 



] p 9/09/86 PROGRAM TO FIT BILL'S 2 POINT RATIONAL EQUATION. CHECKOUT 
] p VERIFICATION READINGS AND THE VALUES AFTER THE RESLOPE AS WELL . 



p MODEL A = P1+02XC+03XCXC 
'MODEL A = S1+P2XC+S3XCXC' 
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:5] 

:e] 

7] 

:8] 

:q] 

:io] 

:ii] 

:i2] 

: 13 ] 

;i4] 
:i5] 
:is] 
; 17 ] 

:i8] 

; 1 9 ] 
! 20 ] 
! 21 ] 
; 2 2 ] 
; 2 3 ] 

!24] 
: 2 5 ] 
:26] 
!27] 
[28] 
; 2 9 ] 
:3o: 
:3i] 
[32] 
; 3 3 ] 
:34] 
: 3 5 ] 
:36] 
:37] 
!3 8] 
!3 9] 
!40] 
!41] 
!42] 
!43] 
!44! 
!45] 
'46] 
!47] 
:48] 
:49] 

; so] 
; 5 1] 

!52] 

: 53 ] 

:54] 
!55] 
!56] 
! 5 7] 
!58] 
!59] 
!60] 
!61] 
!62] 
!63 ] 
!64] 
!65] 
!66] 
: 6 7 ] 
; 6 8 ] 
!69] 
!70] 

171] 

172] 

173] 



P<- 1 0.1 0.25 0.1 0.4 0.25 0.3 0.2 0.2 
E<-20 + 5x x gxj-e-l 
T5<r 8 3 5 6 1 pO 
'FOR TUBE 401 ' 

LOOP:A + + / 3 2p(i ,PPM401 [l+t6;E[J]+i5] )*2 
□<-£<-( 6)4 )1 ( 2 0 5 0 100 xff [I] )• . * 0 1 2 
A ORIGINAL COEFFICIENTS FOR THE TWO CURVES. 

,PPM4 0 1 [ 7 + 1 1 2 ; E [I] + 1 5 ] 

^((-E[2]42x/t:r3] )<100xP[J] ) [OTHER 

O+U-KZ 2] ) + ( (XC2]*2)-4xX[3]xaci] -D [ 2x 1 6 ] ) )*0 . 5 ) * 2 *K [ 3 ] 



<561 
[1] Ul 



0 



C/5<- 8 3 
T5+T5 

C/<-( 6 0*5 61 pL [”l + 2x^6] ) , 8 3 $ 6 1 pO 
' FOR ELEMENT : ' , (EL [I ; ] ) , ' .PRIOR TO RE SLOPE ' 

□<-C/<-( 3 0 ® 6 1 p 10 40 70 xtf [ (p [I] * 5 ) -4] ) ,U 
C<-+/ ( $ ,PPM40 1 [20 + 1 2 ;E[T]+i5])*2 
C^Ct(E+.x(50xP[J] )* 0 1 2 ) 

p GIVES MULTIPLICATIVE CORRECTION FACTORS FOR THE TWO LINES 
L<-<c,PPM401[(22 + i6 ) , 3 4+ i 6 ; E [I] + i 5 ] 

DOQ:DD<rUCZ2l*2)-nxKZ31x(KZll-DZ2xin ) 

^(LL <0 )/ERRORQ 

PCJJ]^Cx((-7f[2] )+DD*0 .5)t2x^[3] 

LQ:+( (II+II+l )56 )/DOQ 
'AFTER THE RESLOPE' 

U+( 6 0 <5 6 1 pPC _ l + 2xi6] ) , 8 3 <5 6 1 po 
D<-$ PPM40 1 [ ( 2 8 + 1 6 )jE[J]+i5] 

II<-1 

KK+ 2 4 4 6 

D01:DD<rOCZ21*2 )-4xtf[3]x(K[l]-£[KK[2T]] ) 

•*(LL < 0 ) / ERROR1 

0[iT]<-Cx((-P[2] )+LL*0 ,5)*2xX[3] 

LI '. + ((11+11 + 1 )S4 )/L 01 

A 0 ^Cx((-K[ 2 ] ) + ( | (X[2]*2)-4x/f[3]x(X[l]-P[2 4 4 6 ] ))*0 
C/<-( 3 0 + C/), [ 1 ] ( ( (6 05 41 pP[l 3 3 5] ), 8 3541 pOK [ 1 ] 
□<-C/<-( 3 0 <5 10 1 p 10 40 70 20 50 50 100 10 40 70 xP[(£T “ 
+((I+I+l)Z9)/LOOP 
I<-1 

LOOT-.A++/ 3 2 p U.PPM401 [41 + 1 6 , E[I]+x 5] 

0+K<-(ts>A) 1(20 50 100 xHZIl )°.* 0 12 
L<-<i> ,PPM401[47 + il2;E[J]+i5] 

■*((-K[2]+2xtf[3] )< 100 xff[J] ) /OTHER1 



) + 2 



1 , Li J 

:i]*5 



.5)*2x/f[3] 
3 0 +£/) 

)-4] ),U 



0 CURV 
i ] ))^0 
:i]-p[ 



p ORIGINAL COEFFICIENTS FOR THE TWO 
-*( ( (XC2]*2)-4xX[3]x(^[i]-D[2x l 6 - 

0<-((-K[2] )+((P[2]*2)-4xP[3]x(p[ 

U 5-e- 8 3 5 6 1 pO 
T5+T5 , [1] i/5 

f/^(6 0 5 6 1 pD[ l + 2xi3]), 8 3 5 6 1 pO 
' FOR ELEMENT : ' , (EL [I ; ] ) , ' , PRIOR TO RE SLOPE ' 
OC/<-(3 0 $ 6 1 p 10 40 70 xfl[(E[I]*5 )-4] ),i/ 
-*-((I<-I+l)£9 )/ LOOT 
■>0 

OPPEPl J £/5<- 8 3 <£.6 1 pO 
C/<- (6 0 5 6 1 pD[ l + 2xi6] ), 8 3 5 6 1 pO 
P5-S-P5 , [1] C/5 

' POP ELEMENT '. 1 , (EL CJ; ]), ' .PRIOR TO RESLOPE ' 
□<-C/<- ( 3 0 5 6 1 p 10 40 70 xP[(P[J] + 5 )-4] ),C/ 

+ ((I<rI+lj£9)/LOOT 
EPP0P1 : 0 [JT] <-8 8 8 8 
+L1 

ERRORQ:OZII]+ 8888 
+LQ 

pEPPOPI :0 [J/]<-88 88 
po-Ll 

p ERRORQ : O C III <-8 8 8 8 
p +LQ 

OTHER :J+I 

U+ (6 0 5 6 1 pL C1+2 xi6]), 8 3 5 6 1 p 0 
C/5<- 8 3 5 6 1 p 0 



C/PKES. 

J/OPPEPl 

2x i 6] ) )*0 . 5 )* 2 xP[ 3 ] 
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:74] 
!7 5 ] 
! 76 ] 

: 7 7 ] 

! 78 ] 

: 7 9 ] 
:so] 
:si] 
! 82 ] 
:s3] 

! 84 ] 

! 85 ] 

:i] 
:2] 
!3 ] 
! 4 ] 
!5] 
! 6 ] 
! 7 ] 

:s] 

!9] 

:io] 

:ii~ 

:i2] 

:i3] 

:i4] 

:i5] 

:i6 _ 

:i7] 

:i8 _ 

:i9 

:20] 

.211 

: 22i 

\23~ 

! 24 ] 

:25] 

;26] 

!27 ] 

!28 

:29 

:3o 

.31 

’.32 

!33 



TS+Tl i, Cl] US 

' FOR ELEMENT : 1 , (EL [J ; ] ) , ' , PRIOR TO RE SLOPE ' 
□<■(/«■( 3 0?s61pl0 40 70 xf/[(E[ ! /]f5)-4]),i/ 

D+$ .PPAJ401C (22 + x 6 ), 3 4 + i6 -,EZJl + \ 5] 

' AFTER THE RE SLOPE ' 

U+(6 0 $ 6 1 pZ? C l+2*x6] ) , 8 3 $ 6 1 pO 
D+to ± PPMn 0 1 C ( 2 8 + 1 6 ) ; E [ J 1 + 1 5 ] 

U+{ 3 0 *£/) , [1] ( ( (6 0 $ 4 1 pj9Cl 3 3 5] ) 



OtM3 0 ? l'O 1 p 10 40 
■*( (J*I+1 )^9 ) /LOOP 
I+-1 
+LOOT 

V SUMOFSQ 

AA+Q1 ,T1 ,T2 ,T3 ,TV,T5 



8 3 541 pO), Cl] 3 
70 20 50 50 100 10'40 70 *H C (EIJ1 t5 )■ 



0 + £/) 
•4] ),U 



E <-& x 1 9 x J-«-i 

LOOP J S<- 12 6 p A A C£’Ci"]+x6;] , Cl] A A C ( 54 +Tx6 ) + x 6 ; ] 
R* 12 1 p$ ,i?l CfiCJ] 6 ; ] 

Ofi2<- 9 2 vR, C2] S 
AW<-7 p 0 
SSUM+7 p 0 
7«-7p0 
I/V^l 

NEST:KO+ 1 
DOiSS+SLINiKOl 
RR+RZINxl 
+(SS> 0 ) /COUNT 
+((KO+KO+l)Z6)/DO 
COUNT :NNZK01 +NN ZKOl+1 
SSUMZK01+SSUMLK01+(SS-RR)*2 
+ ((KO<rKO+l)Z6)/DO 
JK+ 1 

RATIO :^ZNNZJK1-Q )/GO 
VZJK1 +SSUM ZJK1 *NN ZJK1 
+ ( ( JK+JK+1 )£6 ) /RATIO 
+ UIN+IN+ 1 )Z12 ) /NEST 
GO : VZJK1 -<-888 
■*( (JK+JK+1 )<,& ) /RATIO 
+ ( ( IN+IN+l )<>12 ) /NEST 



U+TT+ELZI ;] , ' 
Zl+UU+'SIZE ' , 

Q+VV* ' RATIO ' , 

O ' flAWT ' 

! - 

+ ( (J>I+1 )£9 ) /LOOP 



' IQ 2 * l4SStfM 
10 J 5 J.47V7V 
10 2$ 14-K 

, 10 0 $Mi + 1 + V 
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APPENDIX D 

SUMMARY OF RANK 



DATAAL PPM201 



FE 


175.00 


190.49 


65.59 


119.33 


972.52 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


29.17 


31.75 


10.93 


19.89 


162.09 


888.00 


RANK 




3 


1 


2 


4 


5 


AG 


20.00 


22.11 


33.74 


22.85 


37.23 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


3.33 


3.69 


5.62 


3.81 


6.20 


888.00 


RANK 


* 


1 


3 


2 


4 


5 


AL 


36.00 


49.29 


47.87 


47.76 


47.82 


48.75 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


6.00 


8.22 


7.98 


7.96 


7.97 


8.12 


RANK 


* 


5 


3 


1 


2 


4 


CR 


10.00 


10.04 


9.46 


9.58 


21.59 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


1.67 


1.67 


1.58 


1.60 


3.60 


888.00 


RANK 




3 


1 


2 


4 


5 


CU 


28.00 


29.44 


33.00 


31.86 


42.77 


33.90 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


4.67 


4.91 


5.50 


5.31 


7.13 


5.65 


RANK 


* 


1 


3 


2 


5 


4 


MG 


32.00 


18.54 


20.96 


21.13 


117.57 


21.28 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


5.33 


3.09 


3.49 


3.52 


19.59 


3.55 


RANK 




1 


2 


3 


5 


4 


NI 


19.00 


25.62 


370.96 


14.67 


89.70 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


3.17 


4.27 


61.83 


2.44 


14.95 


888.00 


RANK 




2 


4 


1 


3 


5 


SI 


283.00 


.00 


383.29 


193.42 


.00 


.00 


SIZE 


6.00 


.00 


6.00 


6.00 


.00 


.00 


RATIO 


47.17 


888.00 


63.88 


32.24 


888.00 


888.00 


RANK 




3 


2 


1 


4 


5 


TI 


22.00 


21.31 


167.05 


18.08 


21.07 


19.45 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


3.67 


3.55 


27.84 


3.01 


3.51 


3.24 


RANK 




4 


5 


1 


3 


2 





DATAAL 


PPM202 










FE 


87.00 


90.85 


69.46 


62.50 


387.26 


47.50 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


14.50 


15.14 


11.58 


10.42 


64.54 


7.92 


RANK 




4 


3 


2 


5 


1 


AG 


5.00 


4.67 


3.31 


4.80 


35.16 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


.83 


.78 


.55 


.80 


5.86 


888.00 


RANK 




2 


1 


3 


4 


5 


AL 


11.00 


14.13 


16.47 


16.22 


20.50 


18.22 



48 



SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


1.83 


2.36 


2.75 


2.70 


3.42 


3.04 


RANK 


* 


1 


3 


2 


5 


4 


CR 


2.00 


1.33 


2.03 


1.78 


4.37 


2.28 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


.33 


.22 


.34 


.30 


.73 


.38 


RANK 




1 


3 


2 


5 


4 


CU 


42.00 


38.05 


32.38 


34.15 


59.89 


25.46 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


7.00 


6.34 


5.40 


5.69 


9.98 


4.24 


RANK 




4 


2 


3 


5 


1 


MG 


52.00 


46.76 


47.40 


47.26 


89.28 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


8.67 


7.79 


7.90 


7.88 


14.88 


888.00 


RANK 




1 


3 


2 


4 


5 


NI 


7.00 


4.04 


6486.57 


9.90 


40.81 


15.54 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


1.17 


.67 


1081.10 


1.65 


6.80 


2.59 


RANK 




1 


5 


2 


4 


3 


SI 


304.00 


302.19 


293.80 


289.35 


296.75 


298.05 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


50.67 


50.37 


48.97 


48.23 


49.46 


49.67 


RANK 




5 


2 


1 


3 


4 


TI 


8.00 


8.79 


10.58 


9.53 


9.41 


9.47 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


1.33 


1.46 


1.76 


1.59 


1.57 


1.58 


RANK 


* 


1 


5 


4 


2 


3 





DATAAL 


PPM203 










FE 


230.00 


240.99 


96.31 


179.67 


396.15 


123.58 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


38.33 


40.17 


16.05 


29.95 


66.02 


20.60 


RANK 




4 


1 


3 


5 


2 


AG 


2.00 


1.59 


.53 


1.14 


17.87 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


.33 


.26 


.09 


.19 


2.98 


888.00 


RANK 




3 


1 


2 


4 


5 


AL 


52.00 


45.26 


44.56 


46.30 


58.88 


41.04 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


8.67 


7.54 


7.43 


7.72 


9.81 


6.84 


RANK 




3 


2 


4 


5 


1 


CR 


1.00 


.53 


.48 


.25 


4.27 


.29 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


.17 


.09 


.08 


.04 


.71 


.05 


RANK 




4 


3 


1 


5 


2 


CU 


245.00 


344.41 


386.51 


367.05 


354 . 90 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


40.83 


57.40 


64.42 


61.18 


59.15 


888.00 


RANK 


* 


1 


4 


3 


2 


5 


MG 


9.00 


5.04 


7.18 


7.31 


29.65 


12.93 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


1.50 


.84 


1.20 


1.22 


4.94 


2.15 


RANK 




1 


2 


3 


5 


4 


NI 


9.00 


9.17 


120.78 


17.14 


22.59 


19.61 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 
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RATIO 


1.50 


1.53 


20.13 


2.86 


3.77 


3.27 


RANK 


* 


1 


5 


2 


4 


3 


SI 


305.00 


757.69 


552.89 


1687.68 


753.81 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


50.83 


126.28 


92.15 


281.28 


125.63 


888.00 


RANK 


* 


3 


1 


4 


2 


5 


TI 


13.00 


13 . 54 


14.30 


14.49 


14.32 


14.36 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


2.17 


2.26 


2.38 


2.42 


2.39 


2.39 


RANK 


* 


1 


2 


5 


3 


4 





DATAAL 


PPM206 










FE 


725.00 


823.52 


739.28 


807.14 


979.47 


798.93 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


60.42 


68.63 


61.61 


67.26 


81.62 


66.58 


RANK 


* 


4 


1 


3 


5 


2 


AG 


3.00 


19.31 


17.25 


17.95 


47.74 


15.83 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.25 


1.61 


1.44 


1.50 


3.98 


2.64 


RANK 


* 


3 


1 


2 


5 


4 


AL 


317.00 


111.67 


265.15 


104.02 


6031.69 


.00 


SIZE 


12.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


26.42 


18.61 


44.19 


17.34 


1005.28 


888.00 


RANK 




2 


3 


1 


5 


4 


CR 


7.00 


19.88 


20.78 


20.50 


33.57 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


.58 


1.66 


1.73 


1.71 


2.80 


888.00 


RANK 


* 


1 


3 


2 


4 


5 


CU 


114.00 


170.55 


169.35 


169.37 


183.95 


160.21 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


9.50 


14.21 


14.11 


14.11 


15.33 


13.35 


RANK 


* 


4 


2 


3 


5 


1 


MG 


77.00 


163.79 


172.05 


163.16 


179.81 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


6.42 


13.65 


14.34 


13.60 


14.98 


888.00 


RANK 


* 


2 


3 


1 


4 


5 


NI 


3.00 


93.12 


333.76 


93.53 


162.60 


73.45 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.25 


7.76 


27.81 


7.79 


13.55 


12.24 


RANK 


* 


1 


5 


2 


4 


3 


SI 


748.00 


6691.70 


1407.42 


345.57 


1493.15 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


62.33 


557.64 


117.28 


28.80 


124.43 


888.00 


RANK 




4 


2 


1 


3 


5 


TI 


27.00 


87.04 


78.46 


102.10 


89.66 


52.04 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


2.25 


7.25 


6.54 


8.51 


7.47 


8.67 


RANK 


* 


2 


1 


4 


3 


5 



DATAAL PPM207 


FE 


2898.00 


2914.98 


19613.44 


3395.53 


1457.12 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


6.00 


.00 


RATIO 


241.50 


242.92 


1634.45 


282.96 


242.85 


888.00 


RANK 


* 


2 


5 


3 


1 


4 



50 



AG 


14.00 


154.71 


.26 


52.41 


1.33 


.27 


SIZE 


12.00 


12.00 


6.00 


12.00 


6.00 


6.00 


RATIO 


1.17 


12.89 


.04 


4.37 


.22 


.04 


RANK 




5 


1 


4 


3 


2 


AL 


320.00 


.00 


234.39 


.00 


.00 


.00 


SIZE 


11.00 


.00 


12.00 


.00 


.00 


.00 


RATIO 


29.09 


888.00 


19.53 


888.00 


888.00 


888.00 


RANK 




2 


1 


3 


4 


5 


CR 


11.00 


60.80 


32.88 


52.55 


28.17 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


6.00 


.00 


RATIO 


.92 


5.07 


2.74 


4.38 


4.69 


888.00 


RANK 


* 


4 


1 


2 


3 


5 


CU 


29.00 


290.09 


135.19 


943.80 


20.49 


23.40 


SIZE 


12.00 


12.00 


12.00 


12.00 


6.00 


6.00 


RATIO 


2.42 


24.17 


11.27 


78.65 


3.41 


3.90 


RANK 


* 


4 


3 


5 


1 


2 


MG 


89.00 


395.16 


97.60 


794.64 


107.73 


.00 


SIZE 


12.00 


12.00 


6.00 


12.00 


6.00 


.00 


RATIO 


7.42 


32.93 


16.27 


66.22 


17.95 


888.00 


RANK 


* 


3 


1 


4 


2 


5 


NI 


15.00 


226.95 


271.83 


221.93 


201.59 


4.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


1.25 


18.91 


22.65 


18.49 


16.80 


.67 


RANK 




4 


5 


3 


2 


1 


SI 


443.00 


382.24 


703.85 


528.94 


405.50 


.00 


SIZE 


12.00 


6.00 


12.00 


6.00 


6.00 


.00 


RATIO 


36.92 


63.71 


58.65 


88.16 


67.58 


888.00 


RANK 


* 


2 


1 


4 


3 


5 


TI 


67.00 


230.65 


67.83 


358.10 


171.87 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


5.58 


19.22 


5.65 


29.84 


14.32 


888.00 


RANK 


* 


3 


1 


4 


2 


5 





DATAAL 


PPM302 










FE 


179.00 


295.18 


284.07 


279.17 


1040.62 


327.22 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


14.92 


24.60 


23.67 


23.26 


86.72 


27.27 


RANK 


* 


3 


2 


1 


5 


4 


AG 


12.00 


13.72 


11.42 


12.52 


27.56 


.24 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


1.00 


1.14 


.95 


1.04 


2.30 


.04 


RANK 




4 


2 


3 


5 


1 


AL 


191.00 


27.67 


184.49 


37.31 


.00 


.00 


SIZE 


10.00 


6.00 


6.00 


6.00 


.00 


.00 


RATIO 


19.10 


4.61 


30.75 


6.22 


888.00 


888.00 


RANK 




1 


3 


2 


4 


5 


CR 


6.00 


6.31 


6.63 


6.36 


13.35 


5.10 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


.50 


.53 


.55 


.53 


1.11 


.43 


RANK 




2 


4 


3 


5 


1 


CU 


32.00 


56.40 


64.10 


63.24 


85.96 


65.02 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


2.67 


4.70 


5.34 


5.27 


7.16 


5.42 


RANK 


* 


1 


3 


2 


5 


4 


MG 


85.00 


118.39 


118.44 


114.24 


191.42 


.00 



51 



SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


7.08 


9.87 


9.87 


9.52 


15.95 


888 . 00 


RANK 


* 


2 


3 


1 


9 


5 


NI 


20.00 


21.67 


97.51 


19.37 


100.50 


10.37 


SIZE 


12.00 


12.00 


6.00 


12.00 


12.00 


6.00 


RATIO 


1.67 


1.81 


16.25 


1.20 


8.37 


1.73 


RANK 




3 


5 


1 


9 


2 


SI 


65.00 


97.88 


92.17 


93.20 


92.57 


87.92 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


5 . 42 


8.16 


7.68 


7.77 


7.71 


7.29 


RANK 


* 


5 


2 


9 


3 


1 


TI 


39.00 


9-8.27 


99.77 


51.86 


99.15 


98.79 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


3.25 


9.02 


9.15 


9.32 


9.10 


9.06 


RANK 


* 


1 


9 


5 


3 


2 





DATAAL 


PPM303 










FE 


108.00 


107.97 


111.92 


112.21 


397.53 


123.26 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


18.00 


18.00 


18.65 


18.70 


57.92 


20.59 


RANK 




1 


2 


3 


5 


9 


AG 


.00 


.21 


.36 


.32 


7.19 


.89 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


.00 


.03 


.06 


.05 


1.20 


.19 


RANK 


* 


1 


3 


2 


5 


9 


AL 


6.00 


11.92 


19.39 


11.16 


23.00 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


1.00 


1.99 


2.90 


1.86 


3.83 


888.00 


RANK 


* 


2 


3 


1 


9 


5 


CR 


1.00 


2.29 


3.35 


2.85 


12.25 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


.17 


.38 


.56 


.97 


2.09 


888.00 


RANK 


* 


1 


3 


2 


9 


5 


CU 


90.00 


39.39 


99.98 


53.20 


93.62 


92.99 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


6.67 


6.56 


7.91 


8.87 


7.27 


7.17 


RANK 




1 


9 


5 


3 


2 


MG 


32.00 


29.80 


91.08 


95.75 


39.19 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


5.33 


9.13 


6.85 


7.63 


5.69 


888.00 


RANK 




1 


3 


9 


2 


5 


NI 


6.00 


6.50 


.00 


15.91 


69.70 


.00 


SIZE 


6.00 


6.00 


.00 


6.00 


6.00 


.00 


RATIO 


1.00 


1.08 


888.00 


2.57 


10.78 


888.00 


RANK 


* 


1 


9 


2 


3 


5 


SI 


53.00 


69.80 


505.69 


66.91 


60.38 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


8.83 


10.80 


89.27 


11.15 


10.06 


888.00 


RANK 


* 


2 


9 


3 


1 


5 


TI 


11.00 


11.60 


6.72 


12.59 


12.18 


12.21 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


1.83 


1.93 


1.12 


2.09 


2.03 


2.09 


RANK 




2 


1 


5 


3 


9 



DATAAL1 PPM304 



52 



FE 


64.00 


2180.05 


3250.60 


7185.16 


50624.67 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


10.67 


363.34 


541.77 


1197.53 


8437.45 


888.00 


RANK 


* 


1 


2 


4 


5 


3 


AG 


10.00 


10.03 


9.97 


9.97 


13.44 


8.62 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


1.67 


1.67 


1.66 


1.66 


2.24 


1.44 


RANK 




4 


2 


3 


5 


1 


AL 


72.00 


88.29 


88.09 


87.99 


103.49 


90.89 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


12.00 


14.71 


14.68 


14.66 


17.25 


15.15 


RANK 


* 


3 


2 


1 


5 


4 


CR 


4.00 


3.97 


3.07 


3.32 


13.88 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


.67 


.66 


.51 


.55 


2.31 


888.00 


RANK 




3 


1 


2 


4 


5 


CU 


23.00 


25.91 


24.46 


24.46 


24.99 


26.71 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


3.83 


4.32 


4.08 


4.08 


4.16 


4.45 


RANK 


* 


4 


1 


2 


3 


5 


MG 


78.00 


82.36 


.00 


92.85 


104.65 


.00 


SIZE 


6.00 


6.00 


.00 


6.00 


6.00 


.00 


RATIO 


13.00 


13.73 


888.00 


15.48 


17.44 


888.00 


RANK 


* 


1 


4 


2 


3 


5 


NI 


3.00 


2.32 


22.76 


5.26 


197.67 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


.50 


.39 


3.79 


.88 


32.95 


888.00 


RANK 




1 


3 


2 


4 


5 


SI 


17.00 


16.86 


22.34 


12.79 


16.26 


16.24 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


2.83 


2.81 


3.72 


2.13 


2.71 


2.71 


RANK 




4 


5 


1 


3 


2 


TI 


24.00 


26.31 


27.28 


27.70 


27.08 


26.45 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


4.00 


4.39 


4.55 


4.62 


4.51 


4.41 


RANK 


* 


1 


4 


5 


3 


2 





DATAAL 


PPM305 










FE 


402.00 


460.36 


421.86 


427.68 


3416.65 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


67.00 


76.73 


70.31 


71.28 


569.44 


888 . 00 


RANK 


* 


3 


1 


2 


4 


5 


AG 


2.00 


3.41 


3.68 


3.63 


6.13 


4.19 


SIZE 


5.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


.40 


.57 


.61 


.60 


1.02 


.70 


RANK 


* 


1 


3 


2 


5 


4 


AL 


35.00 


44.39 


21.72 


71.60 


.00 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


.00 


.00 


RATIO 


5.83 


7.40 


3.62 


11.93 


888.00 


888.00 


RANK 




2 


1 


3 


4 


5 


CR 


2.00 


1.87 


2.05 


1.91 


3.59 


2.19 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


.33 


.31 


.34 


.32 


.60 


.36 


RANK 




1 


3 


2 


5 


4 



53 



cu 


165.00 


179.86 


176.91 


179.08 


212.23 


176.33 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


27.50 


29.98 


29.49 


29.85 


35.37 


29.39 


RANK 


* 


4 


2 


3 


5 


1 


MG 


118.00 


118.06 


119.86 


119.84 


467.95 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


19.67 


19.68 


19.98 


19.97 


77.99 


888.00 


RANK 


* 


1 


3 


2 


4 


5 


NI 


35.00 


32.50 


27.84 


21.81 


81.19 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


5.83 


5.42 


4.64 


3.64 


13.53 


888.00 


RANK 




3 


2 


1 


4 


5 


SI 


19.00 


14.18 


17.28 


16.45 


19.92 


16.95 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


RATIO 


3.17 


2.36 


2.88 


2.74 


3.32 


2.83 


RANK 




1 


4 


2 


5 


3 


TI 


38.00 


51.08 


1687.75 


61.65 


51.55 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


6.33 


8.51 


281.29 


10.28 


8.59 


888.00 


RANK 


* 


1 


4 


3 


2 


5 





DATAAL 


PPM401 










FE 


100.00 


128.03 


121.08 


109.09 


3119.08 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


8.33 


10.67 


10.09 


9.09 


259.92 


888.00 


RANK 


* 


3 


2 


1 


4 


5 


AG 


2.00 


.85 


1.70 


1.48 


12.34 


.67 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.17 


.07 


.14 


.12 


1.03 


.11 


RANK 




1 


4 


3 


5 


2 


AL 


268.00 


317.88 


300.21 


297.57 


429.29 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


22.33 


26.49 


25.02 


24.80 


35.77 


888.00 


RANK 


* 


3 


2 


1 


4 


5 


CR 


16.00 


17.29 


17.16 


17.10 


39.32 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


1.33 


1.44 


1.43 


1.42 


3.28 


888.00 


RANK 


* 


3 


2 


1 


4 


5 


CU 


106.00 


102.05 


100.06 


99.74 


100.09 


107.97 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


8.83 


8.50 


8.34 


8.31 


8.34 


9.00 


RANK 




4 


2 


1 


3 


5 


MG 


91.00 


170.50 


156.21 


148.12 


179.10 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


7.58 


14.21 


13.02 


12.34 


14.92 


888.00 


RANK 


* 


3 


2 


1 


4 


5 


NI 


38.00 


34.87 


26.02 


31.70 


192.41 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


3.17 


2.91 


2.17 


2.64 


16.03 


888.00 


RANK 




3 


1 


2 


4 


5 


SI 


162.00 


149.09 


141.29 


166.29 


171.89 


103.04 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


13.50 


12.42 


11.77 


13.86 


14.32 


17.17 


RANK 




2 


1 


3 


4 


5 


TI 


135.00 


142.77 


364.12 


141.42 


134.41 


141.47 



54 



SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


11.25 


11.90 


30.34 


11.79 


11.20 


11.79 


RANK 




4 


5 


2 


1 


3 





DATAAL 


PPM404 










FE 


30.00 


38.60 


92.84 


73.23 


1614.32 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


5.00 


6.43 


15.47 


12.21 


269.05 


888.00 


RANK 


* 


1 


3 


2 


4 


5 


AG 


9.00 


7.92 


6.83 


8.49 


98.37 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


1.50 


1.32 


1.14 


1.42 


16.39 


888.00 


RANK 




2 


1 


3 


4 


5 


AL 


231.00 


427.14 


240.35 


327.70 


966.05 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


38.50 


71.19 


40.06 


54.62 


161.01 


888.00 


RANK 


* 


3 


1 


2 


4 


5 


CR 


1.00 


.97 


2.22 


1.10 


10.77 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


.17 


.16 


.37 


.18 


1.80 


888.00 


RANK 




1 


3 


2 


4 


5 


CU 


113.00 


116.47 


125.29 


127.72 


238.24 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


18.83 


19.41 


20.88 


21.29 


39.71 


888.00 


RANK 


* 


1 


2 


3 


4 


5 


MG 


47.00 


44.59 


45.87 


46.35 


46.32 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


7.83 


7.43 


7.64 


7.73 


7.72 


888.00 


RANK 




1 


2 


4 


3 


5 


NI 


11.00 


8.71 


5.84 


6.41 


145 . 33 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


1.83 


1.45 


.97 


1.07 


24.22 


888.00 


RANK 




3 


1 


2 


4 


5 


SI 


36.00 


53.14 


42.51 


69.86 


820.92 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


6.00 


8.86 


7.08 


11.64 


136.82 


888.00 


RANK 


* 


2 


1 


3 


4 


5 


TI 


32.00 


26.54 


25.41 


26.86 


25.17 


.00 


SIZE 


6.00 


6.00 


6.00 


6.00 


6.00 


.00 


RATIO 


5.33 


4.42 


4.23 


4.48 


4.20 


888.00 


RANK 




3 


2 


4 


1 


5 





DATAAL 


PPM405 










FE 


239.00 


223.84 


175.23 


193.22 


4758.95 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


19.92 


18.65 


14.60 


16.10 


396.58 


888.00 


RANK 




3 


1 


2 


4 


5 


AG 


24.00 


19.51 


18.68 


19.07 


23.69 


14.31 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


2.00 


1.63 


1.56 


1.59 


1.97 


1.19 


RANK 




4 


2 


3 


5 


1 


AL 


120.00 


136.21 


133.36 


135.06 


138.79 


133.60 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


10.00 


11.35 


11.11 


11.25 


11.57 


11.13 


RANK 


* 


4 


1 


3 


5 


2 



55 



CR 


3.00 


6.97 


4.91 


6.10 


14.73 


3.51 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.25 


.58 


.41 


.51 


1.23 


.59 


RANK 


x 


3 


1 


2 


5 


4 


CU 


47.00 


41.02 


47.51 


46.43 


48.57 


49.51 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


3.92 


3.42 


3.96 


3.87 


4.05 


4.13 


RANK 




1 


3 


2 


4 


5 


MG 


266.00 


254.16 


220.37 


194.30 


213.89 


74.36 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


22.17 


21.18 


18.36 


16.19 


17.82 


12.39 


RANK 




5 


4 


2 


3 


1 


NI 


26.00 


32.24 


8.76 


20.00 


167.12 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


2.17 


2.69 


.73 


1.67 


13.93 


888.00 


RANK 




3 


1 


2 


4 


5 


SI 


65.00 


71.11 


71.03 


71.21 


85.72 


60.57 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


5.42 


5.93 


5.92 


5.93 


7.14 


5.05 


RANK 




3 


2 


4 


5 


1 


TI 


40.00 


44.50 


45.33 


43.82 


44.29 


25.18 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


3.33 


3.71 


3.78 


3.65 


3.69 


4.20 


RANK 


x 


3 


4 


1 


2 


5 





DATAAL 


PPM408 










FE 


469.00 


547.44 


363.16 


415.18 


17220.25 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


39.08 


45.62 


30.26 


34.60 


1435.02 


888.00 


RANK 




3 


1 


2 


5 


4 


AG 


33.00 


28.08 


26.34 


26.91 


37.04 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


2.75 


2.34 


2.19 


2.24 


3.09 


888.00 


RANK 




3 


1 


2 


4 


5 


AL 


74.00 


63.13 


58.35 


61.43 


149.71 


10.28 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


6.17 


5.26 


4.86 


5.12 


12.48 


1.71 


RANK 




4 


2 


3 


5 


1 


CR 


2.00 


3.48 


3.37 


3.36 


96.86 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


.17 


.29 


.28 


.28 


8.07 


888.00 


RANK 


x 


3 


2 


1 


4 


5 


CU 


32.00 


35.41 


41.01 


42.99 


43.21 


45.57 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


2.67 


2.95 


3.42 


3.58 


3.60 


3.80 


RANK 


X 


1 


2 


3 


4 


5 


MG 


91.00 


83.41 


95.59 


94.42 


280.57 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


7.58 


6.95 


7.97 


7.87 


23.38 


888.00 


RANK 




1 


3 


2 


4 


5 


NI 


30.00 


28.84 


8.25 


16.98 


478.03 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


2.50 


2.40 


.69 


1.41 


39.84 


888.00 


RANK 




3 


1 


2 


4 


5 



56 



SI 


132.00 


329.10 


306.22 


305.16 


316.09 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


11.00 


27.42 


25.52 


25.43 


26.34 


888.00 


RANK 


x 


4 


2 


1 


3 


5 


TI 


35.00 


32.77 


27.75 


29.39 


27.98 


28.11 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


2.92 


2.73 


2.31 


2.45 


2.33 


2.34 


RANK 




5 


1 


4 


2 


3 





DATAAL 


PPM409 










FE 


915.00 


896.53 


673.04 


788.52 


2643.86 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


76.25 


74.71 


56.09 


65.71 


220.32 


888.00 


RANK 




3 


1 


2 


4 


5 


AG 


41.00 


147.67 


5523.28 


292.54 


90.84 


7.25 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


3.42 


12.31 


460.27 


24.38 


7.57 


1.21 


RANK 




3 


5 


4 


2 


1 


AL 


277.00 


283.60 


265.18 


264.74 


255.95 


265.68 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


23.08 


23.63 


22.10 


22.06 


21.33 


22.14 


RANK 




5 


3 


2 


1 


4 


CR 


3.00 


4.47 


5.43 


5.10 


26.16 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


.25 


.37 


.45 


.42 


2.18 


888.00 


RANK 


x 


1 


3 


2 


4 


5 


CU 


134.00 


146.62 


157.18 


154.99 


218.25 


153.89 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


11.17 


12.22 


13.10 


12.92 


18.19 


12.82 


RANK 


X 


1 


4 


3 


5 


2 


MG 


73.00 


77.80 


78.61 


76.53 


95.99 


66.93 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


6.08 


6.48 


6.55 


6.38 


8.00 


11.15 


RANK 


X 


2 


3 


1 


4 


5 


NI 


3.00 


2.85 


73.03 


6.58 


97.48 


7.02 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.25 


.24 


6.09 


.55 


8.12 


1.17 


RANK 




1 


4 


2 


5 


3 


SI 


525.00 


.00 


216.03 


.00 


5087.06 


.00 


SIZE 


12.00 


.00 


12.00 


.00 


12.00 


.00 


RATIO 


43.75 


888.00 


18.00 


888.00 


423.92 


888.00 


RANK 




3 


1 


4 


2 


5 


TI 


29.00 


25.84 


26.30 


23.93 


26.04 


21.23 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


2.42 


2.15 


2.19 


1.99 


2.17 


3.54 


RANK 




2 


4 


1 


3 


5 





DATAAL 


PPM412 










FE 


41.00 


50.60 


142.39 


76.18 


1098.40 


85.99 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


3.42 


4.22 


11.87 


6.35 


91.53 


14.33 


RANK 


X 


1 


3 


2 


5 


4 


AG 


6.00 


5.50 


3.93 


4.56 


11.75 


1.68 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


.50 


.46 


.33 


.38 


. 98 


.14 



57 



RANK 



2 



3 



5 



1 



4 



AL 


85.00 


241.77 


357.39 


280.80 


53.34 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


6.00 


.00 


RATIO 


7.08 


20.15 


29.78 


23.40 


8.89 


888.00 


RANK 


* 


2 


4 


3 


1 


5 


CR 


3.00 


2.53 


1.94 


2.47 


15.28 


.25 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.25 


.21 


.16 


.21 


1.27 


.04 


RANK 




4 


2 


3 


5 


1 


CU 


33.00 


35.82 


47.30 


44.36 


51.29 


47.38 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


2.75 


2.99 


3.94 


3.70 


4.27 


3.95 


RANK 


* 


1 


3 


2 


5 


4 


MG 


53.00 


42.21 


43.66 


44.11 


163.60 


18.44 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


4.42 


3.52 


3.64 


3.68 


13.63 


3.07 


RANK 




2 


3 


4 


5 


1 


NI 


23.00 


16.25 


809.13 


19.34 


96.26 


22.57 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


1.92 


1.35 


67.43 


1.61 


8.02 


3.76 


RANK 




1 


5 


2 


4 


3 


SI 


880.00 


1541.92 


.00 


.00 


3279.27 


.00 


SIZE 


12.00 


12.00 


.00 


.00 


12.00 


.00 


RATIO 


73.33 


128.49 


888.00 


888.00 


273.27 


888.00 


RANK 


* 


1 


3 


4 


2 


5 


TI 


21.00 


21.31 


20.92 


20.66 


21.39 


20.88 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


1.75 


1.78 


1.74 


1.72 


1.78 


1.74 


RANK 




4 


3 


1 


5 


2 





DATAAL 


PPM413 










FE 


320.00 


298.68 


2303.91 


231.51 


986.67 


147.91 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


26.67 


24.89 


191.99 


19.29 


82.22 


12.33 


RANK 




3 


5 


2 


4 


1 


AG 


17.00 


17.38 


15.35 


16.32 


21.01 


14.15 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


1.42 


1.45 


1.28 


1.36 


1.75 


1.18 


RANK 




4 


2 


3 


5 


1 


AL 


124.00 


267.30 


138.11 


14.47 


16.94 


16.80 


SIZE 


12.00 


12.00 


12.00 


6.00 


6.00 


6.00 


RATIO 


10.33 


22.27 


11.51 


2.41 


2.82 


2.80 


RANK 




5 


4 


1 


3 


2 


CR 


7.00 


4.48 


20.26 


4.47 


15.25 


1.18 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.58 


.37 


1.69 


.37 


1.27 


.20 


RANK 




3 


5 


2 


4 


1 


CU 


81.00 


65.38 


68.28 


61.92 


99.51 


57.06 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


6.75 


5.45 


5.69 


5.16 


8.29 


4.75 


RANK 




3 


4 


2 


5 


1 


MG 


72.00 


66.45 


66.15 


65.62 


110.91 


55.46 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


6.00 


5.54 


5.51 


5.47 


9.24 


9.24 


RANK 




3 


2 


1 


4 


5 



58 



NI 


16.00 


16.21 


306.24 


16,84 


87.17 


14.39 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


1.33 


1.35 


25.52 


1.40 


7.26 


2.40 


RANK 


* 


1 


5 


2 


4 


3 


SI 


548.00 


3253.60 


363.43 


37162.64 


1317.74 


.00 


SIZE 


12.00 


12.00 


6.00 


12.00 


12.00 


.00 


RATIO 


45.67 


271.13 


60.57 


3096.89 


109.81 


888.00 


RANK 


* 


3 


1 


5 


2 


4 


TI 


25.00 


30.57 


26.31 


31.13 


29.47 


29.45 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


2.08 


2.55 


2.19 


2.59 


2.46 


2.45 


RANK 


* 


4 


1 


5 


3 


2 





DATAAL 


PPM415 










FE 


95.00 


98.18 


97.72 


76.81 


768.92 


110.22 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


7.92 


8.18 


8.14 


6.40 


64.08 


9.18 


RANK 




3 


2 


1 


5 


4 


AG 


8.00 


7.57 


7.88 


7.93 


8.29 


7.53 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


.67 


.63 


.66 


.66 


.69 


.63 


RANK 




2 


3 


4 


5 


1 


AL 


67.00 


88.27 


109.28 


88.64 


113.79 


85.55 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


5.58 


7.36 


9.11 


7.39 


9.48 


7.13 


RANK 


* 


2 


4 


3 


5 


1 


CR 


2.00 


1.43 


.93 


.86 


5.70 


.51 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


.17 


.12 


.08 


.07 


.48 


.04 


RANK 




4 


3 


2 


5 


1 


CU 


23.00 


25.87 


27.27 


28.05 


65.71 


32.73 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


1.92 


2.16 


2.27 


2.34 


5.48 


2.73 


RANK 


* 


1 


2 


3 


5 


4 


MG 


31.00 


24.13 


29.92 


30.08 


124.74 


12.40 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


2.58 


2.01 


2.49 


2.51 


10.39 


2.07 


RANK 




1 


3 


4 


5 


2 


NI 


13.00 


12.73 


423.45 


16.17 


39.13 


19.71 


SIZE 


12.00 


12.00 


6.00 


12.00 


12.00 


12.00 


RATIO 


1.08 


1.06 


70.58 


1.35 


3.26 


1.64 


RANK 




1 


5 


2 


4 


3 


SI 


288.00 


197.67 


344.15 


222.31 


174.00 


.00 


SIZE 


11.00 


6.00 


12.00 


6.00 


6.00 


.00 


RATIO 


26.18 


32.94 


28.68 


37.05 


29.00 


888.00 


RANK 


* 


3 


1 


4 


2 


5 


TI 


20.00 


14.81 


9.61 


14.86 


14.62 


.00 


SIZE 


12.00 


12.00 


6.00 


12.00 


12.00 


.00 


RATIO 


1.67 


1.23 


1.60 


1.24 


1.22 


888.00 


RANK 




2 


4 


3 


1 


5 





DATAAL 


PPM416 










FE 


372.00 


1112.62 


906.69 


971 . 94 


1775.91 


93.76 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 



59 



RATIO 


51.00 


92.72 


75.56 


80.99 


147.99 


15.63 


RANK 




4 


2 


3 


5 


1 


AG 


11.00 


23.19 


22.85 


22.93 


18.28 


6.51 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.92 


1.93 


1.90 


1.91 


1.52 


1.09 


RANK 


* 


5 


3 


4 


2 


1 


AL 


81.00 


114.20 


106.24 


108.97 


117.53 


110.48 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


6.75 


9.52 


8.85 


9.08 


9.79 


9.21 


RANK 


* 


4 


1 


2 


5 


3 


CR 


5.00 


19.94 


17.17 


18.37 


18.20 


.89 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.42 


1.66 


1.43 


1.53 


1.52 


.15 


RANK 




5 


2 


4 


3 


1 


CU 


26.00 


141.62 


145 . 05 


140.81 


127.41 


138.40 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


2.17 


11.80 


12.09 


11.73 


10.62 


11.53 


RANK 


* 


4 


5 


3 


1 


2 


MG 


29.00 


106.64 


115.62 


112.50 


121.80 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


2.42 


8.89 


9.64 


9.37 


10.15 


888.00 


RANK 


* 


1 


3 


2 


4 


5 


NI 


15.00 


96.20 


153.68 


90.07 


118.51 


83.23 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


1.25 


8.02 


12.81 


7.51 


9.88 


13.87 


RANK 


* 


2 


4 


1 


3 


5 


SI 


148.00 


928.14 


1098.13 


1078.96 


425.30 


40.81 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


12.33 


77.35 


91.51 


89.91 


35.44 


6.80 


RANK 




3 


5 


4 


2 


1 


TI 


35.00 


53.25 


14.04 


59.28 


57.97 


.00 


SIZE 


12.00 


12.00 


6.00 


12.00 


12.00 


.00 


RATIO 


2.92 


4.44 


2.34 


4.94 


4.83 


888.00 


RANK 




2 


1 


4 


3 


5 





DATAAL 


PPM501 










FE 


180.00 


78.30 


179.09 


131.65 


7817.24 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


15.00 


6.52 


14.92 


10.97 


651.44 


888.00 


RANK 




1 


3 


2 


4 


5 


AG 


10.00 


9.19 


6.30 


7.81 


31.30 


1.95 


SIZE 


11.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.91 


.77 


.53 


.65 


2.61 


.32 


RANK 




4 


2 


3 


5 


1 


AL 


104.00 


139.48 


60.48 


45.26 


46.08 


43.34 


SIZE 


10.00 


12.00 


12.00 


6.00 


6.00 


6.00 


RATIO 


10.40 


11.62 


5.04 


7.54 


7.68 


7.22 


RANK 




5 


1 


3 


4 


2 


CR 


3.00 


5.68 


23.18 


6.15 


31.26 


.00 


SIZE 


11.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


.27 


.47 


1.93 


.51 


2.61 


888.00 


RANK 


* 


1 


3 


2 


4 


5 


CU 


16.00 


11.63 


9.58 


10.06 


19.84 


9.96 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


1.33 


.97 


.80 


.84 


1.65 


.83 



60 



RANK 



1 



3 



5 



2 



4 



MG 


183.00 


187.93 


190.79 


190.92 


1091.13 


.00 


SIZE 


11.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


16.64 


15.66 


15.90 


15.91 


90.93 


888.00 


RANK 




1 


2 


3 


4 


5 


NI 


21.00 


19.30 


141.79 


24.83 


185.73 


.00 


SIZE 


11.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


1.91 


1.61 


11.82 


2.07 


15.48 


888.00 


RANK 




1 


3 


2 


4 


5 


SI 


408.00 


898.40 


499.81 


232.03 


537.23 


564.70 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


34.00 


74.87 


41.65 


19.34 


44.77 


47.06 


RANK 




5 


2 


1 


3 


4 


TI 


23.00 


19.37 


18.96 


18.74 


18.88 


8.07 


SIZE 


11.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


2.09 


1.61 


1.58 


1.56 


1.57 


1.34 


RANK 




5 


4 


2 


3 


1 





DATAAL 


PPM502 










FE 


310.00 


299.60 


267.45 


281.80 


402.95 


260.74 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


25.83 


24.97 


22.29 


23.48 


33.58 


21.73 


RANK 




4 


2 


3 


5 


1 


AG 


4.00 


3.14 


3.51 


4.14 


10.61 


4.13 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


.33 


.26 


.29 


.34 


.88 


.34 


RANK 




1 


2 


4 


5 


3 


AL 


438.00 


488.21 


453.86 


454 . 93 


4541.83 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


36.50 


40.68 


37.82 


37.91 


378.49 


888.00 


RANK 


* 


3 


1 


2 


4 


5 


CR 


6.00 


7.42 


6.21 


7.13 


21.56 


1.12 


SIZE 


11.00 


12.00 


12.00 


12.00 


12.00 


6.00 


RATIO 


.55 


.62 


.52 


.59 


1.80 


.19 


RANK 




4 


2 


3 


5 


1 


CU 


36.00 


40.26 


41.72 


44.04 


48.01 


51.05 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


12.00 


RATIO 


3.00 


3.35 


3.48 


3.67 


4.00 


4.25 


RANK 


* 


1 


2 


3 


4 


5 


MG 


182.00 


169.20 


191.42 


189.20 


2219.15 


.00 


SIZE 


12.00 


12.00 


12.00 


12.00 


12.00 


.00 


RATIO 


15.17 


14.10 


15.95 


15.77 


184.93 


888.00 


RANK 
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